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The  main  objectives  of  the  program  are  to  study  crack  tip  stresses  and  crack 
tip  deformation,  and  to  use  the  findings  of  these  studies  to  investigate  the 
fractures  of  ductile  and  tough  materials. 

Crack  tip  displacement  and  strain  fields  were  measured  with  the  moire  method, 
and  the  crack  tip  stress,  displacement,  and  strain  fields  were  calculated  with 
the  finite  elements  method.  The  measured  displacement  and  strain  fields  agree 
well  with  the  finite  element  calculations.  The  results  were  used  to  develop  the 
direct  correspondence  method  of  evaluating  fracture  toughness  of  ductile  and 
tough  materials  and  were  also  used  to  analyze  deformation  mode  of  fatigue  crack 
growth. 

Griffith^  formulated  his  well  known  energy  criterion  for  brittle  fractures. 
At  the  fracture  initiation  of  a  cracked  brittle  solid  in  the  condition  of  fixed 
grip,  the  release  rate  of  the  stored  strain  energy  equals  or  exceeds  the  dissi¬ 
pative  surface  energy  rate.  3Uc/9a  =  3U$/3a  =  2y,  where  and  U  are  strain 
energy  and  surface  energy;  y,  the  surface  energy  per  unit  area;  and  a,  crack 

length,  y  is  constant  for  a  given  material.  Assuming  a  constant  dissipative 

(2)  (3) 

rate  of  plastic  energy  r,  Irwin'  and  Orowan'  '  extended  the  energy  criterion 
to  metallic  solids,  where  plastic  deformation  takes  place  at  crack  tips. 

The  crack  tip  elastic  stresses,  strains,  and  displacements  are  characterized 

*  } 

by  the  stress  intensity  factor,  K.  Li-.  has  shown  that  for  the  case  of  small 
scale  yielding,  SSY,  K  characterizes  crack  up  stresses,  strains,  and  displace¬ 
ments  even  within  a  crack  tip  plastic  zone  in  a  metallic  specimen.  Liu  has  em¬ 
phasized  that  the  capability  of  K  to  characterize  crack  tip  stresses,  strains  and 
displacements  forms  the  fundamental  basis  of  the  linear  elastic  fracture  mechanics 
rather  than  the  global  energy  balance. 


More  recently,  Hutchinson^ ,  and  Rice  and  Rosengren^  derived  the  charac- 
eristic  crack  tip  stress,  strain  and  displacement  fields  in  power  law  strain 
hardening  materials.  The  crack  tip  stress,  strain  and  displacement  fields  can 
be  characterized  by  J,  which  is  a  contour  independent  integral.  J  is  also  the 
rate  of  potential  energy  change  during  the  cracking  process  in  a  non-linear 
elastic  solid.  J  has  been  widely  used  to  study  non-linear  fracture  mechanics. 

However,  the  characteristic  crack  tip  fields  are  not  universally  valid.  Figure  1 
shows  the  slip  line  fields  of  a  double-edge-notched  specimen  and  a  center-notched 
specimen  under  fully  plastic  tensile  loads.  Because  of  the  difference  in  in-plane 
constraints,  the  principle  tensile  stress  equal  to  (2  +  ir)k  in  one  case  and  2k  in 
the  other;  k  is  the  yield  stress  in  terms  of  shear.  The  characteristic  flow  fields 
and  stress  fields  of  these  two  cases  are  grossly  different.  In  other  words,  for 
different  types  of  specimen  geometry,  the  characteristic  fields  might  be  different. 
The  same  J-value  might  mean  different  crack  tip  stresses,  strains,  and  displacements. 
As  a  result,  one  would  expect  that  specimens  of  different  geometric  types  would 
fail  at  different  J-values.  Indeed,  Hancock  and  Cowling^  tested  specimens  of 
five  different  geometries.  Fig.  2.  The  lowest  and  the  highest  values  of  crack 
tip  opening  displacements  at  fracture,  5C,  differ  by  a  factor  of  10.  6  is  re¬ 

lated  to  J  but  the  relation  differs  for  different  flow  fields.  Taking  the  dif- 

(8) 

ference  of  the  flow  fields  into  consideration  Hancock  and  Cowling'  found  the 
corresponding  extreme  values  of  Jc  to  be  570  kNm'^  and  147  kNm~\  a  difference 
of  a  factor  of  four. 

Compact  tension  specimens  and  three  point  bending  specimens  were  extensively 
used  to  measure  fracture  toughnesses  of  a  number  of  tough  and  ductile  materials. 

Yet  according  to  the  results  shown  in  Fig.  2,  the  measured  Jjj.  values  would  be 
overvalued  in  comparison  with  the  deeply  doublc-edge-cracked  specimens. 
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In  view  of  the  divergent  results  of  recent  findings,  a  sound  and  rational 
method  of  fracture- toughness  evaluation  is  badly  needed.  This  research  pro¬ 
gram  follows  a  much  more  careful  process  of  development .  The  crack  tip  deforma¬ 
tions  were  measured  and  calculated  numerically.  The  measurements  and  the  numeri¬ 
cal  calculations  were  compared  with  the  analytical  results.  Based  on  these 
studies,  the  characteristic  crack  tip  fields  were  deduced.  Then  the  direct 
correspondence  method  of  evaluating  fracture  toughnesses  was  developed. 

Three  characteristic  two-dimensional  crack  models  exist:  plane  strain, 
plane  stress  and  Dugdale  model.  The  characteristic  crack  tip  fields  of  plane 
stress  and  Dugdale  models  were  studied  extensively  in  this  investigation.  The 
crack  tip  stress,  strain  and  displacement  fields  of  a  small  sample  in  large 
scale  yielding  or  general  yielding  correspond  directly  to  the  crack  tip  fields 
of  a  large  sample  in  small  scale  yielding.  Therefore  these  two  samples,  one 
small  and  one  large,  must  have  the  same  K  or  J  value.  This  conclusion  is  sub¬ 
stantiated  by  the  excellent  agreement  between  the  measured  crack  tip  strain 
and  displacement  fields  with  the  finite  element  method  calculations;  the  ex¬ 
cellent  agreement  between  the  finite  element  calculations  and  the  analytical 
results  of  Hutchinson^  and  Rice  and  Rosengren^;  and  the  agreement  of  the 
measured  and  calculated  load-elongation  curves  of  cracked  samples.  The  choice 
of  a  specific  model  and  specific  mechanical  parameter  for  fracture  toughness 
evaluation  is  heavily  dependent  on  the  specimen  thickness.  The  details  are 
summarized  in  Appendix  I:  "Thickness  Effects  on  the  Choice  of  Fracture  Criteria" 

At  times,  the  extensive  experimental  and  theoretical  studies  appear  to  be 
overly  cautious.  However,  in  view  of  the  recent  work  by  Hancock  and  Cowling^, 
our  careful  approach  is  warranted.  The  method  of  direct  correspondence  is  ex¬ 
tremely  promising  for  the  future  development  of  non-linear  fracture  mechanics. 
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Fatigue  cracks  in  ductile  materials  are  often  propagated  by  shear  deforma¬ 
tion  process.  When  a  crack  is  small  and  the  applied  stress  approaches  the  tensile 
yield  stress  of  a  material,  the  crack  tip  plastic  zone  size  will  be  comparable  or 
even  larger  than  the  crack  length.  The  effects  of  large  scale  yielding  will  cer¬ 
tainly  make  the  correlation  between  AK  and  da/dN  invalid.  Crack  tip  shear  defor¬ 
mation  in  both  small  scale  yielding  and  large  scale  yielding  were  calculated. 

Again  the  principle  of  direct  correspondence  is  used  to  establish  the  equivalent 
K-value  of  a  micro-crack  in  large  scale  yielding.  The  results  were  used  to  analyze 
fatigue  crack  growth. 

Materials  are  not  homogeneous.  They  consist  of  weak  and  soft  phases  and  strong 
and  hard  phases.  The  crack  tip  shear  decohesion  displacement  in  a  soft  phase  was 
calculated.  The  shear  deformation  at  the  crack  tip  in  a  soft  phase  is  constrained 
by  the  neighboring  hard  phase.  The  calculated  shear  deformations  were  used  to  analyze 
fatigue  crack  growth.  The  analyses  on  micro-crack  and  the  heterogenous  two-phase 
materials  are  summarized  in  Appendix  II,  "Crack  Tip  Deformation  and  Fatigue  Crack 
Growth." 

The  list  of  publications  on  the  research  results  of  this  program  is  shown 
in  Appendix  III . 

The  grant  has  supported  three  graduate  students:  Drs.  Wan-liang  Hu,  C.  Y. 

Yang  and  Mr.  Chien-Erh  Hong.  Mr.  Hong  has  yet  to  complete  his  thesis.  Their 
dissertation  and  thesis  titles  and  present  affiliations  are  given  in  Appendix  IV. 
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Fig.  2  Crack  tip  openings  at  fracture  initiation  in  range  of  specimen 
geometries  (Ref.  8). 
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ABSTRACT 


The  stresses  and  strains  in  a  cracked  solid  are  in  a  complicated  3-dimensional 
state.  There  are  three  two-dimensional  limiting  cases:  plane  strain,  plane  stress, 
and  Dugdale's  strip  yielding  case.  The  thickness  of  a  plate  relative  to  crack  tip 
plastic  zone  size  determines  which  one  of  these  three  limiting  cases  is  predomi¬ 
nant.  The  characteristics  of  the  plane  stress  and  the  Dugdale  strip  yielding  crack 
tip  fields  were  investigated  with  the  moire  method  together  with  the  finite  ele¬ 
ment  calculations.  The  state  of  the  crack  tip  stress  and  strain  fields  approaches 
that  of  plane  strain  for  very  thick  plates.  For  a  very  thin  plate  with  a  very 
large  crack  tip  plastic  zone,  the  Dugdale's  strip  yielding  model  is  applicable, 
and  the  crack  opening  displacement,  crack  tip  opening  displacement,  and  thickness 
contraction  are  related  to  K  or  J  and  they  can  be  used  for  fracture  toughness 
measurements.  When  a  plate  thickness  is  in  between  these  two  extremes,  there 
exists  a  characteristic  plane  stress  crack  tip  field.  The  correspondence  of  the 
plane  stress  crack  tip  fields  in  small  scale  yielding  and  in  general  yielding 
was  established,  and  the  J-value  of  a  sample  in  general  yielding  can  be  obtained 
from  the  established  correspondence.  The  value  of  the  applied  J-field  of  a  small 
sample  in  general  yielding  can  be  measured  in  terms  of  the  applied  stress  and  the 
specimen  elongation 


Keywords:  Non-linear  fracture  mechanics,  fracture  toughness,  small  scale 
yielding,  general  yielding,  COD,  CTOD,  K,  J,  and  thickness 
contraction. 
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exist.  Still  further  away  is  Region  IV,  where  the  effects  of  specimen  geometry  and 
type  loading  dominate  the  deformation  characteristics .  Note  that,  in  general,  not 
all  of  these  four  regions  are  present.  For  example.  Region  III  diminishes  with 
Increasing  load.  For  a  very  thick  single-edge-cracked,  SEN,  specimen,  only  Region 
I  and  Region  IV  may  exist. 

The  nature  of  the  plane  strain  stiffening  effect  in  Region  I  has  been  discussed 
previously^ It  is  the  stresses  and  strains  in  the  immediate  vicinity  of  the 
crack  tip  within  the  stiffened  zone  that  cause  fracture  initiation.  The  crack  tip 
stiffened  zone  is  embedded  within  the  plane  stress  crack  tip  field.  It  is  reasonable 
that  the  stresses  and  strains  in  the  stiffened  zone  should  be  characterized  by  the 
intensity  of  the  plane  stress  crack  tip  field.  In  other  words,  the  same  intensity 
of  the  plane  stress  crack  tip  field  infers  the  same  stresses  and  strains  in  the 
stiffened  zone,  even  though  the  exact  values  of  the  stresses  and  strains  in  the 
stiffened  zone  are  unknown.  Therefore,  the  intensity  of  the  plane  stress  crack  tip 
field  also  characterizes  the  condition  for  crack  initiation.  The  degree  of  the 
crack  tip  strengthening  is  affected  by  plate  thickness,  therefore  this  conclusion 
Is  applicable  only  to  samples  of  same  thickness.  We  will  focus  our  attention 
in  Regions  II  and  III. 

Crack  tip  stresses,  strains  and  displacements  in  elastic  solids  can  be  expressed 

fol 

In  terms  of  the  stress  intensity  factor  K.  In  a  like  manner,  Hutchinson1  J  and  Rice 
rgl 

and  RosengrenL  J  have  shown  that  crack  tip  stresses,  strains,  and  displacements  can 
be  expressed  in  terms  of  J. 
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II.  CHARACTERISTIC  PLANE  STRESS  CRACK  TIP  FIELD  -  SMALL  SCLE  YIELDING  AND  LARGE 

SCLAE  YIELDING  CORRESPONDENCE 

Crack  tip  strains,  z  t  were  measured  by  Gavigan,  Ke  and  Liu^  in  double-edge- 
cracked,  DEN,  specimens  made  of  three  batches  of  2024-0  aluminum  alloy.  102  mm  and 
203  mm  (4  or  8  inches)  wide  specimens  of  three  different  thicknesses  were  tested. 

The  meeasurements  are  shown  in  Fig.  1.  The  solid  lines  are  plane-stress  finite 
element  calculations.  Away  from  the  crack  tip,  in  the  region  r  >  t,  the  measure¬ 
ments  and  the  calculations  agree  very  well.  Close  to  the  crack  tip,  r  <  t,  the 
measured  strains  are  less  than  the  calculated  values.  These  strains  were  measured 
on  the  specimen  surface;  however  the  surface  was  not  free  to  deform,  because  the 
surface  was  coupled  to  the  interior  of  the  specimen,  where  "plane  strain"  con¬ 
straint  existed.  The  constraint  "stiffened"  the  crack  tip  area  and  it  reduced 
strains  on  the  surface.  The  constraint  is  a  localized  three-dimensional  phenomenon. 
Beyond  the  crack  tip  stiffened  zone,  the  measurements  clearly  indicate  the  existence 
of  a  characteristic  plane  stress  crack  tip  field.  Figure  2  shows  the  measured  and 
the  calculated  load-displacement,  P-A,  curves.  Again,  the  agreement  of  the  measure¬ 
ments  and  the  calculations  is  very  good. 

Figure  3  schematically  shows  various  regions  of  stress  and  strain  fields. 

Region  I  extends  from  the  crack  tip  to  a  distance  approximately  equal  to  half  of 
the  specimen  thickness.  In  this  region,  the  plane  strain  stiffening  effect  exists, 
and  the  three  dimensional  effect  is  clearly  shown.  In  Region  II,  the  plane  stress 
condition  prevails.  The  stress  and  strain  fields  approach  the  analytical  results 
of  Hutchinson^  and  Rice  and  Rosengren^*  Region  III  lies  further  away  from 
the  crack  tip  where  the  material  deforms  elastically  and  the  stresses  and  strains 
vary  with  r'0'5.  The  measurements  shown  in  Fig.  1  were  all  made  at  general  yield¬ 
ing.  However  if  a  specimen  is  large  enough,  the  characteristic  elastic  field  may 

♦This  will  be  discussed  in  detail  later. 
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exist.  Still  further  away  is  Region  IV,  where  the  effects  of  specimen  geometry  and 
type  loading  dominate  the  deformation  characteristics .  Note  that,  in  general,  not 
all  of  these  four  regions  are  present.  For  example.  Region  III  diminishes  with 
increasing  load.  For  a  very  thick  single-edge-cracked,  SEN,  specimen,  only  Region 
I  and  Region  IV  may  exist. 

The  nature  of  the  plane  strain  stiffening  effect  in  Region  I  has  been  discussed 
previously^-1  It  is  the  stresses  and  strains  in  the  immediate  vicinity  of  the 
crack  tip  within  the  stiffened  zone  that  cause  fracture  initiation.  The  crack  tip 
stiffened  zone  is  embedded  within  the  plane  stress  crack  tip  field.  It  is  reasonable 
that  the  stresses  and  strains  in  the  stiffened  zone  should  be  characterized  by  the 
intensity  of  the  plane  stress  crack  tip  field.  In  other  words,  the  same  intensity 
of  the  plane  stress  crack  tip  field  infers  the  same  stresses  and  strains  in  the 
stiffened  zone,  even  though  the  exact  values  of  the  stresses  and  strains  in  the 
stiffened  zone  are  unknown.  Therefore,  the  intensity  of  the  plane  stress  crack  tip 
field  also  characterizes  the  condition  for  crack  initiation.  The  degree  of  the 
crack  tip  strengthening  is  affected  by  plate  thickness,  therefore  this  conclusion 
is  applicable  only  to  samples  of  same  thickness.  We  will  focus  our  attention 
in  Regions  II  and  III. 

Crack  tip  stresses,  strains  and  displacements  in  elastic  solids  can  be  expressed 

in  terms  of  the  stress  intensity  factor  K.  In  a  like  manner,  Hutchinson^  and  Rice 
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and  RosengrenL  have  shown  that  crack  tip  stresses,  strains,  and  displacements  can 
be  expressed  in  terms  of  J. 
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I.  INTRODUCTION 

Griffith^  formulated  the  energy  fracture  criterion  for  brittle  solids.  At 
the  onset  of  fracture  initiation,  the  driving  force  (i.e.  the  rate  of  change  of 
the  elastic  strain  energy  of  a  cracking  solid),  is  equal  to  the  energy  dissipation 
rate  (i.e.  the  surface  energy  of  a  brittle  solid),  which  is  a  constant.  Irwin^ 

m 

and  OrowanL  J  extended  the  Griffith  energy  criterion  to  metallic  solids  by  assuming 
the  plastic  energy  dissipation  rate  for  fracture  initiation  was  a  material  constant. 
However,  it  can  be  shown  that  the  plastic  energy  dissipation  rate  of  a  cracking  in¬ 
finite  plate  in  the  condition  of  plane  strain  is  linearly  proportional  to  crack 
length,  2a.  So  does  the  rate  of  change  of  the  elastic  strain  energy  of  the  plate. 
The  fracture  criterion  of  the  global  energy  balance  leads  to  a  constant  fracture 
stress  for  an  infinite  plate.  The  fracture  stress  is  independent  of  crack  length. 
This  conclusion  contradicts  the  experimental  evidences,  therefore  the  criterion  of 

the  global  energy  balance  for  fracture  initiation  without  the  consideration  of  the 
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detailed  fracture  processes,  must  be  fortuitous1-  ’  . 

Liu'-4,5-^  has  shown  that  the  stress  intensity  factor,  K,  characterizes  crack 
tip  stresses  and  strains  even  within  the  crack  tip  plastic  zone,  rp,  if  the  condi¬ 
tion  of  small  scale  yielding,  S5Y,  is  satisfied.  The  very  existence  of  this  con¬ 
clusion  allows  us  to  determine  experimentally  the  minimum  specimen  size  for  valid 
fracture  toughness  measurements.  But  the  condition  of  SSY  is  a  sufficient,  not 
the  necessary,  condition  for  the  validity  of  the  linear  elastic  fracture  mechanics. 
The  necessary  condition  is  that  K  should  be  able  to  characterize  the  crack  tip 
stress  or  strain  responsible  for  the  defined  fracture  process.  When  the  applied 
stress  or  strain  at  the  fracture  process  zone  reaches  the  critical  value,  frac¬ 
ture  will  occur. 

If  it  is  limited  to  the  realm  of  the  linear  elastic  fracture  mechanics,  with 
the  restriction  to  the  same  fracture  mode  and  the  limitation  to  the  same  thickness, 
the  concept  of  the  K-characterization  of  the  crack  tip  field  and  the  assumption  of 
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the  constant  dissipation  rate  of  plastic  deformation  energy  lead  to  the  same  re¬ 
sults.  However,  when  it  is  extended  beyond  the  realm  of  the  LEFM,  only  the  correct 
analysis  will  lead  to  the  desired  results.  The  concept  of  characterizing 
crack  tip  stress  and  strain  fields  by  the  stress  intensity  factor,  K,  forms  a  sound 
basis  for  the  extension  of  the  LEFM  to  non-linear  fracture  mechanics.  This  paper 
illustrates  its  applications. 

The  characteristic  crack  tip  stresses  and  strains  are  greatly  affected  by 

specimen  thickness.  When  a  small  plastic  zone  at  a  crack  tip  is  imbedded  in  a 

massive  and  thick  plate,  the  constraint  to  thickness  contraction  induces  a  state  of 

high  triaxial  stresses,  and  the  state  of  the  crack  tip  field  approaches  that  of 

plane  strain.  The  maximum  principle  stress,  omflx  is  nearly  2.45  x  oflow,  and  the 

triaxial  state  of  stresses  restrains  plastic  deformation.  On  the  other  hand,  a 

very  large  plastic  zone  in  a  thin  plate  causes  crack  tip  necking.  In  this  instance, 

extensive  plastic  deformation  is  concentrated  within  a  thin  narrow  strip.  Within 

the  strip,  the  plate  thickness  contracts  considerably,  and  the  deformation  approaches 

rsi 

the  Dugdale  strip  yielding  model1  .  When  specimen  thickness  is  intermediate  between 
these  two  extremes,  a  charcteristic  plane  stress  crack  tip  field  exists.  In  this 
case,  the  maximum  principle  stress  is  nearly  equal  to  afriow- 

Fracture  processes  are  controlled  by  the  crack  tip  stresses  and  strains  which 
are  greatly  affected  by  plate  thickness.  Therefore  the  choice  of  the  fracture 
criteria  must  also  be  dependent  on  plate  thickness. 

This  paper  summarizes  and  synthesizes  the  experimental  and  theoretical  work  on 
the  effects  of  thickness  on  crack  tip  stresses,  on  crack  tip  deformation,  and  on 
the  choice  of  the  fracture  criteria  for  non-linear  fracture  mechanics.  Plate  thick¬ 
nesses  less  than  that  of  plane  strain  were  studied.  The  results  are  applicable  to  t 
structures  made  of  high  toughness  materials  which  are  less  thar  2  inches  thick. 
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for  a  material  with  a  power  law  stress-strain  relation 
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where,  aQ  is  the  flow  stress  at  the  strain  eQ .  When  aQ  =  and  eQ  =  ey ,  a  =  1 . 
n  >  1 .  is  effective  stress,  a..  (0,n),  a  (0,n),  e . .  (0,n)  and  u.  (0,n)  are 
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functions  of  9  and  n.  They  define  the  distributions  of  their  corresponding  com¬ 
ponents  in  the  0-direction.  In(n)  is  a  function  of  n. 
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J  is  the  well  known  J-integral.  It  characterizes  crack  tip  stresses,  strains,  and 
displacements  analogous  to  K  for  elastic  solids. 

For  a  non-linear  elastic  solid,  J  can  also  be  considered  as  the  rate  of  the 
change  of  the  potential  energy  with  respect  to  crack  increment.  In  the  case  of  SSY, 


E 
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where  E  =  E  for  plane  stress  and  E  =  E/(l-v  )  for  plane  strain,  v  is  the  Poisson's 
ratio. 

However,  the  characteristic  crack  tip  fields  given  by  Eqn.  (1)  are  not  uni¬ 
versally  valid  in  all  the  cases.  Figure  (4)  shows  the  slip  line  fields  of  a 
double-edge-notched  specimen  and  a  single-edged-notched  specimen  under  fully  plas¬ 
tic  tensile  loads.  Because  of  the  difference  in  in-plane  constraint,  the  principle 
tensile  stress  equal  to  (2  +  tt ) k  in  one  case  and  2k  in  the  other;  k  is  the  yield 
stress  in  terms  of  shear.  The  characteristic  flow  fields  and  stress  fields  of 
these  two  cases  are  grossly  different.  In  other  words,  for  different  types  of 
specimen  geometry,  the  characteristic  fields  might  be  different.  The  same 
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J-value  might  indicate  different  crack  tip  stresses,  strains,  and  displacements. 

As  a  result,  one  would  expect  that  specimens  of  different  geometric  types  would 
fail  at  different  J-values.  Indeed,  Hancock  and  Cowling^^  tested  specimens  of 
five  different  geometries.  The  lowest  and  the  highest  values  of  crack  tip  open¬ 
ing  displacements  at  fracture,  6c  differ  by  a  factor  of  10.  6  is  related  to  J 

but  the  relation  differs  for  different  flow  fields.  Taking  the  difference  of 
the  flow  fields  into  consideration  Hancock  and  Cowling^1^  found  the  correspond¬ 
ing  extreme  values  of  Jc  to  be  570  kNm'^  and  147  kNm'^ ,  a  difference  of  a  factor 
of  four. 

In  order  to  measure  fracture  toughness  as  a  material  property,  a  reference 
state  of  crack  tip  field  must  be  established.  The  most  meaningful  one  is  that  of 
the  small  scale  yielding.  With  this  choice,  the  measured  fracture  toughness  of 
a  small  sample  in  general  yielding  will  be  able  to  predict  the  fracture  strength 
of  a  large  structure  in  small  scale  yielding. 

Plane  stress  calculations  are  also  made  for  the  double-edge-cracked  sample 
and  the  single-edge-cracked  sample  in  SSY,  LSY,  and  GY.  The  characterization 
of  the  crack  tip  fields,  the  correspondence  between  the  crack  tip  fields  in  the 
state  of  the  SSY,  LSY  and  GY,  and  their  relation  to  the  applied  stress  and  the 
Imposed  elongation  are  studied. 

The  crack  tip  stress  and  strain  fields  in  samples  made  of  two  batches  of  alumi¬ 
num  alloy  in  the  state  of  SSY  were  studied  using  plane  stress  finite  analysis.  A 
circular  sample  with  a  crack  along  one  of  the  radial  lines  and  with  boundary  dis¬ 
placements  given  by  the  dominant  terms  of  the  linear  elastic  solution,  i.e.  u^  * 
(K/2y)(r/2ir)^  1^(9),  serves  as  the  reference  state,  see  Fig.  5.  u  is  shear 
modulus.  The  a  -  e  relations  are  prescribed  by  a  =  Ee  for  o  <  Oy  and  o  =  keN 
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for  a  >  a  y.  k  =  ENay^'N^.  ay  =  69  MPa  (10  ksi)  and  N  =  0.22  for  batch  B;  and 
oy  =  50  MPa  (7.26  ksi)  and  N  =  0.315  for  batch  C.  The  finite  element  mesh  is 
shown  in  Fig.  5.  The  details  of  the  calculations  are  given  in  Ref.  12.  The  crack 
line  stresses  and  strains  for  batch  C  aluminum  are  shown  in  the  dimensionless  plots 
in  Fig.  6.  The  results  of  nine  loading  steps  are  shown.  The  value  of  Kj  ranges 
from  1.24  to  4.4  MPa^T  (1.25  to  4  ksi/TF). 


Within  r  ,  the  y-direction  strain  and  stress  can  be  expressed  as 
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where  r^  =  y(Kj/oy)  =  yJ/eyay.  The  values  of  3,  S'  and  y  are  respectively  1.41, 
0.98,  and  0.243  for  batch  C  aluminum,  and  1.35,  0.92,  and  0.281  for  batch  B 
aluminum. 

The  form  of  Eq.  (9)  is  similar  to  that  given  by  Hutchinson^ ,  and  Rice  and 
[9] 

RosengrenL  ,  Eq.  (1).  Combining  Eqs .  (1  and  9),  with  (m+m1)  =  1,  one  obtains 


In(n)/Oyy(°*n)£yy(°»n)  =  BB  7Y 


(10) 


TABLE  I 

N 

n 

6 

8' 

Y 

BB'/Y 

I  ^nVo  (0,n)e  (O.n) 

n  yy  yy 

Batch  B 

0.22 

4.5 

1 .35 

0.92 

0.28 

4.4 

4.4 

Batch  C 

0.32 

3.2 

1  .41 

0.98 

0.24 

5.7 

5.0 

Table  I  lists  the  values  of  (B8'/y)  obtained  from  the  finite  element  calculation 
and  the  values  of  (I eyy)  obtained  from  Hutchinson's  results.  The  values  of 
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cryy  (0,n)  and  e  (0,n)  are  obtained  by  the  linear  interpolation  between 
cryy  (0,3)  =  1.1  and  ayy  (0,13)  =  1.2;  and  cyy  (0.3)  =  0.7  and  eyy  (0.13)  =  0.8; 
the  values  of  In(n)  are  from  the  linear  interpolation  between  I  (3)  =  3.86  and 
I  (5)  =  3.41.  The  values  of  the  finite  element  calculations  and  those  obtained 
from  Hutchinson's  results  are  very  close. 

The  results  for  a  ,  e  ,  a,  and  ep  along  the  radial  lines  of  45°,  60°,  and 

J  V  </  J 

90°  and  away  from  the  crack  line  are  shown  in  Fig.  7  as  solid  curves,  a  is 
effective  stress;  ep  is  effective  plastic  strain.  Similar  relations  can  also  be 
obtained  for  these  stresses  and  strains. 

The  state  of  the  crack  tip  stress  and  strain  fields  in  SSY  serves  as  a 
reference  to  be  compared  with  the  crack  tip  fields  of  specimens  in  LSY  and  GY. 

Double-edge-notched  (DEN)  and  single-edge-notched  (SEN)  specimens  were 

studied.  The  calculated  a  ,  £yy,  a,  and  op  are  shown  in  Figs.  7  and  8.  For 

comparison,  the  results  of  the  SSY  calculations  (solid  curves)  are  also  shown. 

The  values  of  rp  for  the  specimens  in  LSY  and  GY  are  obtained  by  the  linear 

extrapolation  of  effective  stress,  a,  to  yield  stress  Oy  in  a  logarithmic  plot 

of  a  vs  r.  The  solid  lines  of  the  SSY  coincide  with  the  data  of  GY.  The  same 

scaling  factor  rp  is  applicable  to  all  stress  and  strain  components.  Therefore, 

these  stresses  and  strains  can  be  expressed  in  terms  of  a  single  parameter,  r  . 

Since  r  is  related  to  K  or  J  in  SSY,  these  stresses  and  strains  can  also  be 
P 

written  in  terms  of  K  or  J.  In  other  words,  these  stresses  and  strains  in  LSY 
and  GY  are  characterized  by  K  or  J. 

The  fact  that  the  data  of  SSY  and  GY  fall  onto  the  same  line  means  that  a 
crack  tip  field  in  GY  does  correspond  uniquely  to  a  crack  tip  field  in  SSY.  Con¬ 
sequently,  one  will  be  able  to  use  a  small  sample  in  GY  to  reach  the  crack  tip 
field  of  a  large  sample  in  SSY  at  the  equivalent  K-value.  The  fracture  process 


is  controlled  by  the  crack  tip  stresses  and  strains  within  the  stiffened  zone, 
and  the  crack  tip  stiffened  zone  is  imbedded  within  the  plane  stress  zone.  Thus, 
the  stresses  and  strains  within  the  stiffened  zone  will  be  the  same  if  the  intensity 
of  the  surrounding  plane  stress  crack  tip  field  is  the  same.  Or,  to  state  the  con¬ 
clusion  more  simply,  the  K  or  J  of  the  plane  stress  crack  tip  field  characterizes 
the  stresses  and  strains  within  the  stiffened  zone,  even  though  the  stresses  and 
strains  within  the  stiffened  zone  are  unknown.  Therefore,  if  a  small  sample  in 
GY  fractures  at  a  K-value  so  defined,  a  very  large  sample,  made  of  the  same 
material  in  SSY,  will  fracture  at  the  same  K-value,  because  the  crack  tip  stresses 
and  strains  are  identical  in  both  samples.  Since  crack  tip  stresses  and  strains 
in  the  stiffened  zone  is  affected  by  plate  thickness,  this  conclusion  is  valied 
only  for  specimens  of  the  same  thickness. 

The  near  tip  stresses  and  strains  correlate  with  K  or  J.  In  order  to  use  the 

correlation  for  convenient  fracture  toughness  measurements ,  the  near  tip  field  must 

be  related  to  the  applied  stress  ok  and  the  imposed  elongation  A.  Fig.  9  shows 

the  correlation  of  the  product  a ~  at  a  certain  distance  r  ahead  of  a  crack 

tip  with  the  product  o  A  for  the  DEN  sample.  The  product  a  c  w  is  linearly  pro- 

cn  yy  yy 

portional  to  J.  Therefore,  we  have 

K2 

~  -  J  =  13.4  a  A  (batch  C  aluminum) 

E  00 

(ID 

y2 

=  J  =  13.7  a  A  (batch  B  aluminum) 

E  00 

These  relations  can  be  used  to  measure  the  fracture  toughnesses  of  the  DEN  speci¬ 
mens  made  of  these  two  materials. 

In  a  cracked  elastic  solid,  an  area  exists  near  the  crack  tip  bounded  by 
rg,  within  which  the  singular  terms  of  the  characteristic  crack  tip  stress  and 
strain  fields  are  valid  approximations.  For  metallic  specimens  it  can  be  shown 
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that  if  size  of  rg  »  r  ,  K  will  characterize  crack  tip  stresses  and  strains 

f4  s') 

even  w!'  •  .  rpL  ’  J.  The  very  existence. of  this  conclusion  serves  as  the  basis 

for  thf  iiirical  determination  of  the  minimum  specimen  size  for  valid  fracture 

tough;;  measurements.  The  condition  of  r„  »  r  is  known  as  the  condition  of 

e  p 

small  ...ale  yielding.  Wilson  has  found  that  the  size  of  rg  is  quite  small^3^ 
when  t  vy-'-ed  with  other  specimen  dimensions.  However,  rg  is  linearly  proportional 

to  spt  •  . .  size;  so,  in  principle  the  condition  of  SSY  can  always  be  satisfied 

by  ui;. y  large  enough  sample. 

Vv  condition  of  rg  »  rp  could  be  unduly  restrictive  in  terms  of  specimen 

size  .jirements.  T*'  condition  of  r  >>  r  is  a  sufficient  condition.  The 

e  p 

neces:...  condit/bsjs  validity  of  the  linear  elastic  fracture  mechanics  is 

that  k  ni  ld  be  .  ■<  characterize  the  crack  tip  stress  or  strain  component  at 
the  If-  inn  the  defined  fracture  process. 


,  iyire  10  >,hows  the  relation  between  the  applied  stress  and  the  equivalent 
K-va';'  of  a  double-edge-notched  specimen,  both  in  the  linear  elastic  fracture 

mechanic;  region  and  in  the  non-linear  region.  The  linear  elastic  solution  is 

•* 

capab1-  o  characterize  the  crack  tip  field  up  to  a^/Oy  =  0.6,  which  gives  a 

crack  i  ;e  plastic  zone  size  rp(0_g)  -  0.25  inches  in  comparison  with  the  crack 
lengti,  =  0.8  inches.  In  this  case,  the  condition  of  rfi  »  rp  is  obviously 
more  r.  irictive  than  necessary. 

: ;  can  be  concluded  that  the  plane  stress  crack  tip  field  of  a  double-edge- 
notchf  or  single-edge-notched  specimen  in  GY  correlates  well  with  the  crack  tip 
field  SSY ,  and  that  the  crack  tip  field  of  a  small  sample  in  GY  can  be  ex- 
pressc  i  terms  of  K  or  J.  The  correlation  is  substantiated  by  the  experimental 

strain  asurements  and  the  measured  load-elongation  curve  as  well  as  by  its 

agreeni  with  the  analytical  results  of  Hutchinson  and  Rice  and  Rosengren. 


The  ch^-  teristic  plane  stress  crack  tip  field  is  related  to  the  product  of  the 
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applied  stress  and  the  imposed  elongation.  The  relation  can  be  used  to  measure 
fracture  toughnesses  of  small  samples  in  6Y.  The  size  of  the  crack  tip  stiffened 
zone  is  approximately  equal  to  the  plate  thickness.  In  order  to  use  the  charac¬ 
teristic  plane  stress  field,  the  plastic  zone  perhaps  should  be  more  than  twice 
the  plate  thickness;  i.e.,  0.25  (Kc/Oy)  _>  2t.  For  a  tough  material,  with  K£  = 
352  MPa^in  (320  ksi/Tn),  a y  =  552  MPa  (80  ksi),  t  could  be  up  to  two  inches. 

III.  CRACK  TIP  OPENING  DISPLACEMENT  AS  A  FRACTURE  CRITERION 

When  a  very  large  crack  tip  plastic  zone  exists  in  a  thin  plate,  a  strip 
necking  zone  takes  place,  as  observed  by  Dugdale^.  Plastic  deformation  is  con¬ 
centrated  within  the  strip  necking  zone,  and  the  plate  thickness  is  greatly  re¬ 
duced  within  the  strip.  A  strip  necking  zone  is  imbedded  within  a  much  larger 
plastic  zone,  shown  scehmatically  in  Fig.  11,  and  a  strip  necking  zone  is  shown 
clearly  by  the  moire  pattern  in  Fig.  12. 

F141 

Schaeffer  et.  alL  J  have  measured  crack  opening  displacements  as  well  as 
the  relative  "opening  displacements"  between  the  upper  and  the  lower  boundaries 
of  the  strip  necking  zone.  Their  measurements  agree  very  well  with  the  calculated 
values  according  to  the  Dugdale  strip  yielding  model.  The  length  of  the  strip 
necking  zone  has  to  be  several  times  the  thickness  of  the  plate  in  order  to  assure 
a  good  agreement.  The  length  of  the  strip  necking  zone,  relative  to  plate  thick- 
ness,  t,  is  given  by  the  parameter  n  =  (K/Oy)  /t.  For  the  specimens  tested  by 
Schaeffer  et  al ,  n  =  48.  Furthermore,  they  have  found  that  the  opening  displace¬ 
ments  within  the  strip  necking  zone  equal  the  thickness  contraction. 

As  K  increases,  and  a  crack  grows  slowly  toward  a  point  ahead  of  the  crack 
tip,  the  relative  opening  displacement  at  the  point  in  the  strip  necking  zone  in¬ 
creases  while  the  thickness  is  reduced.  Once  the  crack  tip  passes  the  point,  the 
plate  thickness  at  the  point  stops  changing.  After  a  specimen  is  broken,  the 
thickness  at  a  point  along  the  fracture  path  is  the  same  thickness,  when  the  crack 
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tip  first  reached  the  point.  The  thickness  contraction.  At,  at  the  point  equals 

the  crack  tip  opening  displacement,  CTOD,  and  CTOO  is  directly  related  to  the 

Cl  51 

applied  K.  According  to  the  Dugdale  model L  J, 

CTOD  =  K2/Eay  (12) 

for  small  scale  yielding.  Therefore  the  K-values  during  stable  crack  growth  can 
be  obtained  by  plate  thickness  measurements  along  the  fracture  surface  of  a 
broken  specimen.  Indeed,  Liu  and  Kuo^^  have  measured  the  fracture  resistance 
curve  as  well  as  Kc  from  thickness  measurements  in  such  a  manner.  The  thickness 
contraction  fracture  toughness  measurements  can  be  used  economically  for  screening 
tests. 

When  strip  necking  takes  place,  CTOD  is  the  result  of  the  stretch  of  the 
materials  within  the  strip  necking  zone.  In  this  case,  the  crack  tip  opening  dis¬ 
placement  is  a  "tangible"  physical  quantity,  that  can  be  measured  directly. 

Experimenal  observations  showed  that  crack  tip  necking  began  to  form  at  :i 
equal  to  about  18^17^.  For  specimens  in  which  n  is  between  18  and  48,  the  strip 
necking  zone  is  not  large  enough  to  warrant  the  use  of  the  thickness  contraction 
as  a  fracture  toughness  measurement. 

Using  the  moire  method^7’^,  COD  were  measured  in  specimens  made  of  2024-0, 

2024-T3,  and  2024-T351  aluminum  alloys  with  tensile  yield  strengths,  54,  310,  and 
2 

383  MN/m  ,  respectively.  All  specimens  were  101.6  mm  wide  central  crack  panels 
with  a  nominal  crack  length  of  17.8  mm.  The  applied  stresses  were  about  one  half 
the  yield  strengths  for  all  of  the  specimens.  The  density  of  the  moire  grille  is 
528  lines/mm,  and  is  capable  of  measuring  COD  of  about  2.5  x  10‘4  mm  (10*5  in), 
using  the  rotational  mismatch  technique^19^ .  Figure  13  shows  the  measured  COD. 

The  dashed  lines  are  calculated  from  the  Dugdale  model  and  the  solid  lines  from 
the  elastic  model.  Figure  13a  shows  two  specimens  at  the  same  K-value.  Specimen 
1  is  much  thinner  than  specimen  2.  The  thinner  specimen  (n  =  1.75)  agrees  very 
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well  with  the  Dugdale  model,  and  the  thicker  one  (n  =  1.12)  agrees  with  the  elas¬ 
tic  model.  Approaching  the  crack  tip,  the  difference  in  COD  is  more  than  a  factor 
of  two.  It  is  not  likely  that  the  crack  tip  opening  displacement  of  specimen  2 
will  be  equal  to  that  given  by  the  Dugdale  model.  Similar  data  is  shown  in  Fig. 
13b.  It  can  be  concluded  that  when  n  is  close  to  18,  the  measured  COD's  agree 
exceedingly  well  with  the  Dugdale  model.  When  n  is  close  to  one  or  less,  the 
measured  COD's  agree  well  with  the  elastic  calculations.  When  the  COD  measure¬ 
ments  agree  with  the  Dugdale  model,  the  measured  COD  can  be  used  to  infer  a  CTOD 
according  to  the  Dugdale  model  and  the  inferred  CTOD  relates  to  K  by  Eq.  (12). 

The  inferred  CTOD  characterizes  crack  tip  deformation  and  crack  tip  stresses. 
Therefore  it  is  a  suitable  choice  as  a  fracture  criterion.  It  should  be  emphasized 
that  no  strip  necking  zone  was  observed  in  the  thin  specimens  despite  the  fact  that 
the  COD's  in  these  specimens  agreed  well  with  the  Dugdale  model,  which  assumes  a 
strip  yield  zone. 

The  agreement  in  COD  between  the  measurements  and  the  Dugdale  model  was  ob¬ 
served  for  aluminum  specimens  in  small  scale  yielding  or  prior  to  general  yielding. 
Figure  14  shows  the  COD  measurements  in  a  magnesium  specimen  by  Kobayashi  et  al 
All  curves  are  beyond  general  yielding.  The  points  shown  in  the  figure  were 
calculated  in  this  investigation  with  the  Dugdale  model  for  curve  Mg- 1  - 1 Y ,  which 
was  at  the  incipience  of  general  yielding.  The  measurements  agree  well  with  the 
Dugdale  model  even  in  the  general  yielding  condition.  The  agreement  of  the 
measured  and  the  calculated  COD  both  in  SSY  and  in  GY  supports  the  Dugdale  model 
for  fracture  toughness  measurements.  Figure  14  also  shows  that  crack  opening  pro¬ 
files  are  parallel  once  beyond  general  yielding.  When  the  upper  and  the  lower 
crack  surfaces  move  apart  like  two  rigid  surface,  the  area  ahead  of  the  crack  tip 
is  stretched,  and  the  added  stretch  gives  more  opening  displacement  at  the  crack 
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tip.  During  the  process  of  crack  tip  opening  stretching,  the  crack  profile 
maintains  its  original  shape,  corresponding  to  that  at  the  incipience  of  general 
yielding.  Let  6C0D(r)  be  the  difference  of  the  crack  opening  displacement  at  r, 
COD(r),  and  the  crack  tip  opening  displacement  at  the  incipience  of  general  yield¬ 
ing.  Once  beyond  general  yielding,  the  5C0D(r)  remains  unchanged.  The  total 
COD(r)  beyond  general  yielding  is  the  sum  of  <SC0D(r)  and  CTOD,  i.e., 

COD(r)  =  SCOD(r)  +  CTOD 

or  (13) 

CTOD  =  COD(r)  -  6C0D(r) 

6C0D(r)  can  be  calculated  by  the  finite  element  method.  CTOD  can  be  obtained  from 
the  measured  COD(r)  and  the  calculated  6C0D(r).  CTOD  is  related  to  K  by  Eq.  (12). 

Equation  (13)  is  applicable  to  center  cracked  panels.  For  the  compact  ten¬ 
sion  specimens,  an  additional  term  should  be  added  to  account  for  the  rigid  body 
rotation  about  a  "plastic-hinge".  The  rotational  effects  remain  to  be  investigated. 

The  finite  element  method,  based  on  the  Dugdale  model,  can  be  used  to  calcu- 

T21 1 

late  <$C0D(r)  for  plates  of  arbitrary  shape.  Hayes  and  William1-  J  have  used  the 
finite  element  method  to  calculate  COD  and  CTOD  for  the  central  crack,  double¬ 
edge  crack,  single-edge  crack,  and  circumferential  crack  specimens  under  various 
load  levels  up  to  general  yielding.  Basically,  Hayes  and  Williams's  method 
balances  the  stress  intensity  factors  contributed  from  the  remote  applied  stress 
and  the  stresses  in  the  strip  yielding  zone. 

Equation  (12)  has  been  indiscriminantly  used  to  relate  CTOD  to  K  without  the 
consideration  of  the  thickness  effect.  The  equation  is  applicable  when  n  _>  18. 

For  specimens  of  lower  n-values,  Eq.  (12)  needs  to  be  modified  by  a  parameter 
which  is  thickness  dependent. 
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There  exist  several  methods  to  measure  CTOD,  such  as  the  load-line  mouth 
opening  and  plastic  hinge  method,  the  linear  extrapolation  of  a  crack  opening  pro¬ 
file  from  an  optical  photograph,  rubber  infiltration,  etc.  None  of  these  methods 
can  avoid  some  degree  of  arbitrariness,  either  in  the  relation  between  the  far 
field  measurement  and  the  near  tip  deformation  behavior  or  in  the  relation  between 
the  measured  quantity  and  the  inferred  K-value.  Perhaps  some  of  the  earlier  dif¬ 
ficulties  in  applying  CTOD  to  fracture  toughness  measurements  arise  from  such 
arbitrariness.  The  method  and  the  procedure  outlined  in  this  study  will  help 
to  avoid  these  difficulties. 

IV.  SUMMARY  AND  CONCLUSIONS 

1.  The  stresses  and  strains  in  a  cracked  solid  are  in  a  complicated  3-dimensional 
state.  There  are  three  2-dimensional  limiting  cases:  plane  strain,  plane  stress, 
and  the  Dugdale  strip  yielding  case.  The  parameter  n  =  (K/Oy)  /*  determines  which 
one  of  these  three  limiting  cases  is  predominant.  The  characteristic  crack  tip 
fields  of  these  three  cases  differ  vastly. 

2.  Fracture  process  is  controlled  by  the  stresses  and  strains  at  the 
crack  tip.  Because  of  the  difference  in  the  characteristics  of  the  crack 

tip  field,  there  exist  various  regimes  of  fracture  correlations  between  a  small 
laboratory  sample  in  general  yielding  and  a  large  engineering  structure  in  small 
scale  yielding. 

3.  Crack  tip  strip  necking  has  been  observed,  when  n  is  48.  At  n  =  48,  the 
measured  relative  opening  displacements  between  tne  upper  and  the  lower  boundaries 
of  the  strip  necking  zone,  including  that  at  the  crack  tip,  agree  well  with  the 
calculated  values  according  to  the  Dugdale  strip  yielding  model.  Furthermore 

the  relative  opening  displacement  in  the  strip  necking  zone  is  equal  to  the  thick¬ 
ness  contraction.  According  to  the  Dugdale  model, 

K  =  (Eoy  CT0D)1/2  =  (Eay  At)1/2 

where  At  is  thickness  contraction. 
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4.  When  n  is  18,  crack  tip  necking  zone  begins  to  form,  and  the  measured  crack 
opening  displacements  agree  extremely  well  with  the  calculations  based  on  the 
Dugdale  strip  yielding  model.  Once  general  yielding  is  reached,  the  upper  and 
the  lower  crack  surfaces  move  apart  as  two  rigid  surfaces.  The  difference  be¬ 
tween  the  crack  opening  displacement  at  a  distance  r,  COD(r),  and  crack  tip  open¬ 
ing  displacement,  CTOO,  remains  unchanged  once  beyond  general  yielding. 

6C0D(r)  =  COD(r)  -  CTOD 

6C0D(r)  can  be  calculated  by  the  finite  element  method.  CTOD  can  be  obtained 
from  the  calculated  6C0D(r)  and  the  measured  COD(r);  and  CTOD  is  related  to  K. 
Thus  for  n  between  18  and  48,  we  have 

K  =  (Eay  CTOD)1/2 

*  (Eay  [COD(r)  -  SC0D(r)]}1/2  . 

5.  When  rp  is  several  times  larger  than  the  plate  thickness  but  it  is  not  large 
enough  to  form  the  crack  tip  necking  zone,  a  characteristic  plane  stress  crack 
tip  field  exists.  For  double-edge-notched  and  single-edge-notcheJ  samples,  there 
is  a  unique  correspondence  between  the  crack  tip  field  is  «  small  ; vs->l e  in  gen¬ 
eral  yielding  and  the  crack  tip  field  in  a  very  large  sample  in.  small  scale  yield 
ing.  Therefore  the  equivalent  K-value  of  such  a  small  sample  can  be  established. 
Furthermore,  the  near  tip  stresses  and  strains  are  related  to  the  applied  stress 
and  the  imposed  elongation.  Therefore  the  fracture  toughnesses  of  such  small 
samples  can  be  measured  conveniently  for  n  between  10  and  18. 

In  summary,  as  shown  in  Fig.  15,  when  n  is  more  than  48,  thickness  contrac¬ 
tion  can  be  used  to  measure  K;  when  ri  is  between  18  and  48,  the  COD  can  be  used 
to  measure  equivalent  K;  and  when  n  is  between  10  and  18,  the  correspondence  be¬ 
tween  the  crack  tip  fields  in  small  scale  yielding  and  in  general  yielding  can  be 
used  to  measure  the  fracture  toughness  of  a  ductile  and  tough  material.  However, 
the  demarcations  between  the  various  regions  are  tentative.  Additional  studies 
are  needed  to  establish  them  firmly. 
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0.45 
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FIGURE  13.  THICKNESS  EFFECTS  ON  COD. 

Dashed  line  -  Dugdale  model . 
Solid  line  -  Elastic  model. 


FIGURE  14.  COD  IN  MAGNESIUM  [From  A.  S.  Kobayashi  et.  al.  Ref.  20]. 
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FIGURE  15.  THICKNESS  EFFECTS  ON  THE  CHOICE  OF  FRACTURE  CRITERIA. 
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ABSTRACT 

Fatigue  crack  growth  is  caused  by  crack  tip  cyclic  plastic  deformation. 

Both  the  macro-analysis  and  the  moire  strain  measurements  indicate  that  crack 
tip  cyclic  deformation  is  a  function  of  AK,  R-ratio,  and  plate  thickness,  t. 
Therefore,  da/dN  must  be  dependent  on  these  three  parameters. 

The  cyclic  crack  tip  tinzipping  shear  decohesion  process  is  described,  and 
the  unzipping  fatigue  crack  growth  process  is  modelled  by  the  finite  element 
method.  The  calculated  crack  growth  rate  agrees  well  with  the  experimental 
measurements.  The  unzipping  model  is  used  to  study  the  growth  of  micro-cracks 
and  the  fatigue  crack  growth  in  a  ferritic-martensitic  two  phase  steel. 

A  model  of  fatigue  crack  growth  threshold  is  proposed.  The  proposed  model 
agrees  with  the  observed  crack  growth  behavior  in  the  near  threshold  region.  A 
quantitative  analysis  of  fatigue  limit  and  microstructure  is  made.  The  analysis 
provides  a  quantitative  approach  to  optimize  the  microstructure  for  high  fatigue 
strength. 

INTRODUCTION 

Plastic  deformation  in  crystalline  materials  takes  place  in  discrete  paral¬ 
lel  planes  separated  by  crystal  layers  which  are  practically  "deformation  free" 
[l].  This  deformation  process  may  occur  at  a  crack  tip  and  cause  a  crack  to  grow 
under  cyclic  loading. 

The  use  of  a  single  parameter  AK  to  characterize  crack  tip  cyclic  plastic 
deformation  has  been  analyzed  [2,3].  Moire  strain  measurements  [4]  showed  that 
AK  does  characterize  both  the  cyclic  strain  range  and  the  total  accumulated 
strain  in  the  crack  tip  region.  The  capability  of  AK  to  characterize  crack  tip 
deformation  is  the  basis  for  the  use  of  AK  to  correlate  fatigue  crack  growth 
rate.  The  analysis  and  the  moire  strain  measurements  are  reviewed  briefly. 

Alternating  shear  at  two  intersecting  conjugate  slip  planes,  crack  tip 
blunting,  and  shear  decohesion  have  been  proposed  as  the  mechanisms  for  fatigue 
crack  growth  [5, 6, 7, 8].  The  unzipping  model  of  fatigue  crack  growth  based  on 
crack  tip  shear  decohesion  process  was  modelled  by  the  finite  element  method  [9] . 
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The  calculated  crack  growth  rate  compares  favorably  with  crack  growth  rate  and 
strlatlon  spacing  measurements.  The  unzipping  model  was  extended  to  study  the 
growth  of  micro-cracks  and  the  fatigue  crack  growth  in  a  ferritic-martensitic 
two  phase  steel  [10,11]. 

When  the  applied  AK  is  low  enough,  a  crack  stops  propagating.  AK^  is  the 
AK  for  crack  growth  threshold.  Based  on  the  unzipping  shear  decohesion  process, 
a  model  for  AKtjj  is  proposed.  The  proposed  model  agrees  with  the  crack  growth 
behavior  in  the  near  threshold  region. 

A  quantitative  analysis  on  the  fatigue  limit  of  a  two-phase  ferritic-marten¬ 
sitic  steel  is  made.  The  fatigue  limit  of  such  steel  could  be  controlled  by  the 
crack  non-nucleation  or  crack  non-propagation  in  the  ferrite  or  in  the  martensite; 
the  fatigue  limit  could  be  controlled  by  the  constraint  of  the  hard  martensite 
on  the  crack  tip  deformation  in  the  ferrite;  it  can  also  be  controlled  by  the 
cracking  resistance  of  the  martensite.  The  relation  between  the  dominant  con¬ 
trolling  process  and  the  ferrite  domain  size  is  analyzed. 

CRACK  TIP  DEFORMATION  -  AK  CHARACTERIZATION  [2,3] 

During  the  last  twenty  years,  linear  elastic  fracture  mechanics  has  been 
applied  extensively  to  fatigue  crack  growth.  Linear  elastic  fracture  mechanics 
is  based  on  the  characteristic  crack  tip  stress  field  in  a  linear  elastic  solid. 
The  elastic  stresses,  strains,  and  displacements  for  a  Mode  I  crack  are 
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where  r,  0  are  polar  coordinates;  and  a  ,  e  ,  u±  prescribe  the  distributions 
of  the  corresponding  Ojj  ,  £ij ,  and  u^  components.  Kj  is  the  stress  intensity 
factor  while  V  is  Poisson's  ratio. 

These  relations  are  valid  only  in  a  region  very  close  to  a  crack  tip,  and 
we  call  this  region  the  characteristic  elastic  crack  field  zone,  re.  Outside 
this  region,  the  approximations  of  these  relations  deteriorate  rapidly.  Within 
re,  the  stresses  increase  rapidly  as  r  approaches  the  tip.  The  high  stresses 
cause  plastic  deformation  within  a  small  region  close  to  the  tip  in  a  metallic 
specimen.  The  plastic  region  is  defined  by  rp(6). 

Figure  1  shows  two  samples  that  have  different  geometry  and  forms  of  load¬ 
ing.  However,  these  two  samples  are  loaded  to  the  same  K-value.  Without  plas¬ 
tic  deformation,  the  elastic  stresses  within  re's  of  these  two  samples  are  the 
same.  With  plastic  deformation,  if  rp  «  re,  the  stresses  on  the  outer  bound¬ 
aries  of  re's  are  not  perturbed  much  by  the  stress  relaxation  within  rp's,  and 
the  stresses  on  the  boundaries  of  these  two  re's  are  essentially  given  by  the 
elastic  solution. 

Figure  2  shows  the  two  regions  of  re's  as  free  bodies  of  the  same  shape  and 
size.  If  the  condition  of  r_  «  re  is  satisfied,  their  boundary  stresses  must 
be  the  same  since  K^  ■  K2.  Therefore,  the  stresses  and  strains  at  the  same  point 
P(r,8)  in  these  two  free  bodies  also  must  be  the  same.  Furthermore,  if  the  same 
cyclic  loading  (in  terms  of  K) ,  are  applied  to  these  two  samples,  the  fluctua¬ 
tions  of  the  boundary  stresses  on  these  two  free  bodies  are  the  same,  and  the 
cyclic  stresses  and  strains  at  the  same  point  P(r,9)  within  these  two  regions  of 


Fig.  1.  Two  different  specimen  geometry,  two  different 
forms  of  loading,  but  =  K2. 


Fig.  2.  Two  free  bodies  of  same  shape  and  size.  Uith 
same  boundary  stresses. 


re  must  be  the  same  even  within  rp.  If  the  same  fluctuation  of  K  is  applied  to 
these  two  samples  (provide^  the  ratio  rp/a  is  not  excessively  large),  the  cyclic 
stresses  and  strains  experienced  by  these  two  material  points  at  P(r,6)  are  the 
same,  even  if  the  crack  tip  advances.  In  conclusion,  the  fluctuation  of  K  char¬ 
acterizes  the  cyclic  stresses  and  strains  in  the  near  tip  region  of  a  propagat¬ 
ing  crack,  even  within  rp.  This  characterization  is  invariant  to  the  planar  geo¬ 
metry  of  a  specimen  if  tne  condition  of  small  scale  yielding,  SSY  (i.e.  rp «  re) 
is  satisfied. 

When  a  crack  propagates  under  a  constant  AK  (or  under  a  slowly  changing  AK) 
at  a  given  R-ratio,  the  value  of  AK  alone  will  be  able  to  characterize  the  cyclic 
stresses  and  strains  in  the  near  tip  region.  If  the  cyclic  stresses  and  strains 
experienced  by  the  materials  control  the  crack  growth  rate,  then  one  must  con¬ 
clude  that 


da/dN  =  f(AK) 


(2) 


The  stresses  and  strains  are  also  affected  by  specimen  thickness,  t,  and  R-ratio. 
The  functional  relationship  between  da/dN  and  AK  is  affected  by  environment,  cy¬ 
clic  frequency,  V,  and  testing  temperature,  T.  Consequently,  for  a  given  mater¬ 
ial,  one  may  write 


da/dN  *  f (AK,R,t,v,T,  environment)  (3) 

This  relation  is  true  if  r  <<  re  .  re  is  linearly  proportional  to  speci¬ 
men  size.  Therefore,  Equation  13)  can  always  be  made  true,  in  principle,  by 
using  large  enough  samples.  The  SSY  condition  is  the  sufficient  rather  than  the 
necessary  condition.  The  necessary  condition  for  Equation  (3)  is  that  a  sample 
is  large  enough  so  that  K  characterizes  the  specific  stress  or  strain  component 
controlling  the  crack  growth  process. 

The  capability  of  AK  to  characterize  cyclic  strains  in  the  near  tip  region 
was  investigated  by  the  moire  method  and  the  results  are  given  in  the  next 
section. 

CYCLIC  DEFORMATION  IN  CRACK  TIP  REGION  [4] 

The  above  analysis  shows  that  under  a  constant  AK  (or  slowly  changing  AK) 
fatigue  loading,  crack  tip  stresses  and  strains  are  characterized  by  AK.  Be¬ 
sides  AK,  the  stresses  and  strains  are  also  affected  by  specimen  thickness  and 
R-ratio.  Within  the  plastic  zone  in  a  thick  specimen  in  the  state  of  plane 

strain,  the  ratio  between  the  maximum  tensile  stress  and  flow  stress,  (0  /oc,  ) 

max  t  low 

is  approximately  2.45  under  a  static  loading;  and  the  ratio  is  nearly  1  for  the 
state  of  plane  stress.  The  strains  in  these  two  different  states  must  vary 
accordingly. 

The  cyclic  deformations  in  the  crack  tip  regions  at  various  K-levels  were 
measured  with  the  moire  technique  [4] .  Measurements  were  made  on  2024-T351 
aluminum  center  cracked  plate  specimens  with  a  thickness  of  6.4  mm  and  1.28  mm 
(0.25"  and  0.05")  and  R-ratios  of  1/3  and  1/10.  The  static  and  cyclic  yield 
strengths  were  360  and  460  MPa  (52  and  67  ksi)  respectively,  and  the  cyclic 
strain  hardening  exponent  was  0.13. 

As  a  fatigue  crack  propagates  toward  a  material  point  far  ahead  of  the  tip, 
the  cyclic  strain  at  the  point  is,  at  first,  elastic.  As  the  tip  gets  closer  to 
the  point,  the  material  at  the  point  experiences  plastic  deformation  during  the 
loading  half  cycle,  but  the  deformation  is  still  elastic  during  the  unloading 
half  cycle.  The  plastic  strain  accumulates  monotonlcally  cycle  after  cycle  as 
the  tip  propagates  toward  the  point.  As  the  crack  tip  moves  even  closer,  plas¬ 
tic  deformation  takes  place  during  the  unloading  as  well  as  the  loading  half 
cycles,  and  the  plastic  deformation  becomes  cyclic .  In  a  fatigue  loading,  there 
are  monotonic  plastic  zone  r  ^  and  cyclic  plastic  zone,  The  cyclic 

plastic  zone  is  embedded  within  the  monotonic  plastic  zone.  Both  the  accumulated 
total  strain,  emax,  and  the  cyclic  strain  range,  Ae,  along  the  crack  line  were 

measured.  From  these  strain  measurements,  the  monotonic  and  cyclic  plastic 
zones  were  determined  [4]. 

The  moire  patterns  were  obtained  by  the  double  exposure  method.  The  speci¬ 
men  surface  was  coated  with  a  photo-resist  coating.  The  coating  was  first  ex¬ 
posed  to  a  high  density  moire  grille  of  13,000  lpi,  before  the  pre-cracked  sample 
was  loaded.  The  grille  lines  were  parallel  to  the  crack  line  to  measure  the  dis¬ 
placement  u  .  The  second  exposure  was  made  at  K  ^  when  the  crack  increment  was 

long  enough  so  that  a  "steady  state"  of  cyclic  plastic  deformation  at  the  crack 
tip  region  was  realized.  The  moire  fringe  pattern  thus  formed  gave  the  total 

accumulated  strain,  e  .  After  the  picture  of  the  fringe  pattern  was  taken,  the 

max 

pattern  was  immediately  removed,  and  a  new  photo-resist  coating  was  applied  to 


—a sure  cyclic  strain  range,  Ae .  Two  exposures  were  made  in  succession,  one 

each  at  K  .  and  K  .  The  fringe  pattern  gave  Ac. 
min  max 

The  data  of  e  and  Ae  for  the  1/4"  thick  specimens  atR=l/10  are  given  in 

ZD3X 

Figures  3  and  4.  Additional  data  are  given  in  Reference  [4].  The  data  at  a 
given  AK-  or  Kmav-level  in  a  logarithmic  plot  can  be  correlated  very  well  by  a 

straight  line.  For  the  c  data,  the  slopes  of  the  lines  decrease  from  -0.5  at 

max 

low  K  levels  to  -0.8  at  high  K  levels, 
max  max  „ 

Unfortunately,  all  the  Ac  measurements  were  made  outside  the  cyclic  plastic 

zone.  The  moire  pattern  lost  its  definition  when  r  <  0.03".  However,  all  of 

the  Ae  measurements  were  made  within  the  monotonic  plastic  zone.  The  lines  in 

the  log  Ac  vs  log  r  plot  have  a  slope  of  -1/2,  which  agrees  with  the  elastic 

calculation.  If  both  a  and  a  are  taken  into  account, 

xx  yy 

C  ”  (1  -  v)0  /E.  a  (r,0)  ■  K/ /27rr  .  The  measured  Ac  can  be  written  as 

yy  yy  yy  yy 


Ac 


yy  (measured) 


=  gAe 


yy (theoretical) 


(4) 


The  values  of  8  are  close  to  1,  and  they  are  tabulated  in  Table  I.  The  maximum 
deviation  from  the  theoretical  values  is  35%.  For  a  given  combination  of  t  and 
R,  8  is  a  constant. 


t 

1.27  mm  (0.05") 
1.27  mm  (0.05") 
6.35  mm  (0.25") 
6.35  mm  (0.25") 


TABLE  1 

R  6 

1/10 

1/3 

1/10 

1/3 


^measured 

^theoretical 

1.06 

1.18 

1.25 

1.35 


The  measured  Ac  in  the  thin  specimens  is  less  than  that  in  the  thick  speci- 
and  is  caused  by  crack  closure  in  the  thin  specimens  as  will  be  explained  later. 
Ae  is  higher  at  a  higher  R-ratio,  The  higher  Ae  is  likely  caused  by  cyclic 

creep,  which  is  caused  by  the  positive  mean  stress  within  r  ... 

plLJ 


The  cyclic  plastic  zone  size  r  .c.  can  be  obtained  by  extrapolating  the 
measured  Ae  line  to  a  value  of  9.3  8'ld-3,  which  corresponds  to  2oyrQ\  ■ 

938  MPa  (134  ksi)  with  the  Mises  yield  criterion.  The  monotonic  plastic  zone 
can  be  determined  in  a  like  manner.  The  results  are  shown  in  Figures  5  and  6. 
The  dashed  lines  in  these  two  figures  are  the  calculated  lines  based  on  Irwin's 
plastic  zone  size  formula.  The  larger  measured  plastic  zone  sizes  in  comparison 
with  the  calculated  ones  reflect  the  high  measured  strains. 


Recognizing  2oY^  -  yield  stress  range,  at  the  initiation  of  the  cyclic 
plastic  deformation,  we  have, 


u  2a  .  _ ak_ 

^yyOO  e  2<7y(c)  b  e  ^ 


p(C) 


(5) 
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Fig.  7.  Strain  versus  r/r^^j  for  6.35  mm  thick  specimens 
at  R  *  1/10.  Ref.  [4b] 

Similarly,  the  £max  data  are  plotted  against  r/r  in  the  same  figure. 

The  e  data  lie  within  a  narrow  band.  The  empirical  relation 
max 


correlates  well  with  the  data,  t  can  also  be  written  as 

max 

e  “  a  (I— - — — - )  (8) 

max  2tt  2oY(c)  2  Y(C)r 

which  is  similar  to  that  given  by  Hutchinsen  [12],  Rice  and  Rosengren  [13] .  For 
t  -  6.35  mm  and  R  -  1/10,  the  data  in  Figure  7  give  m  -  0.8  and  a  -  0.0045.  Ad¬ 
ditional  data  are  shown  in  Figures  8,  9,  and  10  for  the  other  combination  of  t  and 
R.  e  is  the  highest  for  the  thin  specimens  at  the  high  R-rato  of  1/3,  and  the 

lowest  for  the  thick  specimens  at  the  low  R-ratio  of  1/10.  The  other  two  sets 
of  data  are  in  between.  The  Moire  method  measures  the  deformation  on  the  surface 
of  a  specimen.  The  plastic  deformation  on  the  surface  of  a  thick  specimen  is 
constrained  by  the  interior  of  the  specimen,  where  the  deformation  is  in  the 
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Fig.  10. 


Strain  versus  r/r 

P(C) 

specimens  at  R  *  1/3. 


for  1.27  mm  thick 
Ref.  [4b] 


Fig.  11.  Fatigue  crack  propagation  rate  for  2024-T351 
aluminum  alloy.  Ref.  [4^] 


state  of  plane  strain.  On  the  surface  of  a  thin  specimen,  the  constraint  is  ab¬ 
sent,  and  it  is  "free"  to  deform.  Therefore  e  ,  in  thin  specimens  are  higher  -11- 

QlcUC  S 

than  those  in  thick  specimens.  The  high  measured  e  is  also  reflected  in  the 
i  , -j  £  max 

large  r  ,  .  in  Figure  6. 

pis; 

Both  cyclic  creep  and  crack  closure  will  affect  the  strains.  At  high  R- 
ratio,  the  mean  stress  in  the  plastic  zone  is  positive;  the  positive  mean  stress 
causes  cyclic  creep.  Cyclic  creep  increases  both  At  and  e  at  the  higher  R- 

ratio  of  1/3  as  shown  in  Table  I  and  Figures  8,  9,  and  10  an$  11. 

The  increase  in  emax  is  more  pronounced  in  the  thin  specimens  because  of 
the  lack  of  the  constraint  by  the  interior  plane  strain  zone.  The  material  with¬ 
in  the  plastic  zone  in  thin  specimens  is  stretched  considerably  more  than  that 
in  the  thick  specimens.  During  the  unloading,  the  highly  stretched  zone  will 
cause  crack  closure  and  will  reduce  the  effective  AK.  Therefore,  Ae  in  a  thin 
specimen  is  less  than  that  in  a  thick  specimen  as  shown  in  Table  I. 

Figure  7  shows  that  all  the  Ae  data  fall  on  a  single  line,  therefore,  Ae 
can  be  expressed  in  terms  of  r^^,  which  is  in  turn  related  to  AK.  Hence,  Ae 

can  be  directly  related  to  AK.  The  emax  data  fall  within  a  narrow  band  (±  5%). 

Therefore,  e  can  also  be  expressed  in  terms  of  r  and  AK.  At  a  given 
’  max  v  p(C) 

thickness  and  at  a  given  R-ratio,  AK  can  characterize  both  Ae  and  e 

max 

The  e^x  scatter  bands  for  the  other  combinations  of  t  and  R  are  wider 
(±  15%)  as  shown  in  Figures  8,  9,  and  10.  If  the  lines  for  the  three  high  Up¬ 
values  are  eliminated  from  Figures  8  and  9,  the  scatter  band  is  reduced  con¬ 
siderably.  The  possible  causes  for  data  scatter  include  the  thumb  nail  configu¬ 
ration  of  the  internal  crack  front,  and  the  shift  of  the  crack  to  the  shear  plane 
at  the  high  1^  levels.  The  shift  reduced  plane  strain  constraint  and  caused 
an  increase  in  plastic  strain. 

The  fatigue  crack  growth  data  in  Figure  11  show  that  the  crack  growth  rate 
increases  in  the  same  order  as  the  increase  in  Ae,  that  is  in  the  order  of  R  = 

1/10,  t  =  1.27  mm;  R  =  1/3,  t  =  1.25  mm;  R  =  1/10,  t  «  6.35  mm;  and  R  *  1/3,  t  * 

6.35  mm. 

In  conclusion:  AK  characterizes  both  Ae  and  e  in  a  crack  tip  region  and 

max 

the  cyclic  strains  are  affected  by  both  specimen  thickness  and  R-ratio.  Since 
fatigue  crack  growth  is  closely  related  to  plastic  deformation,  therefore  da/dN 
must  be  related  to  AK,  t,  and  R. 

THE  UNZIPPING  MODEL  OF  FATIGUE  CRACK  GROWTH 

Orowan  [5]  has  proposed  the  alternating  shear  rupture  on  two  sets  of  inter¬ 
secting  slip  planes  as  the  ductile  fracture  mechanism.  Subsequently,  crack  tip 
blunting  by  shear  decohesion  has  been  proposed  as  the  mechanism  for  fatigue  stri- 
ation  formation  and  fatigue  crack  growth  [6,7,8].  Neumann  [8]  has  directly  ob¬ 
served  alternating  shear  rupture,  which  causes  crack  tip  opening  and  crack 
growth  in  copper  single  crystals  under  cyclic  loading. 

Plastic  deformation  in  single  crystals  is  not  continuous  as  has  been  assumed 
in  continuum  mechanics.  Rather,  it  concentrates  In  discrete  parallel  slip  bands 
separated  by  crystal  layers,  which  are  essentially  elastic  as  shown  in  Figure  12. 

The  same  deformation  processes  may  also  occur  at  the  crack  tip  in  a  crystalline 
solid.  The  slip  bands  in  a  single  crystal  coincide  with  the  slip  planes  of  the 
crystal.  The  deformation  in  the  slip  band  of  a  polycrystal  are  the  result  of  a 
number  of  active  slip  systems  of  each  crystal  in  a  band,  and  the  slip  bands  co¬ 
incide  with  the  plane  of  maximum  shear  stresses.  Figure  13  shows  two  sets  of 
parallel  intersecting  slip  bands  at  a  crack  tip.  As  the  applied  stress  increases, 
shear  decohesion  takes  place  alternately  on  slip  bands  1,  2,  3,  and  4,  in 
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Fig.  12.  Development  of  laminar  slip  bands  with  Increasing  defor¬ 
mation,  (a)  and  (b),  and  laminae  at  high  deformation  (c). 
Ref.  [1]. 
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Fig.  13.  The  unzipping  model  for  crack  opening  and  advance. 


succession.  The  elastic  slabs  between  the  neighboring  slip  bands  move  away  from 
the  crack  tip,  one  at  a  time,  like  the  teeth  of  a  zipper  during  the  unzipping 
process.  As  the  slabs  move,  the  crack  tip  opens  up,  and  the  crack  moves  forward. 

Kuo  and  Liu  [9]  have  modelled  the  crack  tip  unzipping  process  using  the 
finite  element  method.  A  cracked  circular  solid  is  loaded  with  boundary  dis¬ 
placement  exitatlon  as  shown  in  Figure  14a.  The  solid  is  loaded  incrementally 
to  a  predetermined  K-level.  The  stresses  and  strains  are  calculated,  and  the 
plastic  zone,  r„,  is  delineated.  At  the  K-level,  a  slip  band  such  as  AB  is 
chosen.  Figure  14b.  Each  of  the  nodal  points  along  AB  is  branched  into  two.  A 
small  load  increment,  6K  is  then  applied.  During  the  increment,  the  branched 
nodal  points  are  allowed  to  slide  freely  along  AB,  but  they  are  not  allowed  to 
move  away  from  AB.  In  the  meantime,  plastic  deformation  is  allowed  to  occur  to 
all  the  elements  during  6K.  The  crack  increment  caused  by  the  unzipping  process, 

6Aa  ,  during  the  increment  6K  is  calculated, 
uz 

According  to  the  dimensional  analysis  and  the  results  of  the  past  studies, 
the  unzipping  crack  increment  must  be 


2  K 

Aa  „  -  C(l-v  )  £- 
uz  Eo,, 


(9) 


where  C  is  an  unknown  proportional  constant;  v  is  the  Poisson's  ratio;  E  is  the 
Young's  modulus ;and  is  the  yield  stress.  Differentiating  both  sides,  one 
obtains 


Fig.  14a.  Finite  element  configu¬ 
ration,  $  =  45°. 


Fig.  14b.  Plastic  zone  and  decohesion 
plane. 
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6Aa  =  C(l-V  )  — — 

uz  Eo„ 


(10) 


The  values  of  K,  6K,  and  £Aa  from  the  finite  element  calculation  enable  us  to 

9  w  UZ 

evaluate  C.  Thus  we  have 


Aa 

uz 


0.04  (1-v2)  f-  =  ^ 
Eay  aY 


(11) 


2  2 

where  J  =  (1  -  v  )K  /E.  Yang  [14]  repeated  the  calculations  for  a  cracked  cir¬ 
cular  solid  and,  in  addition,  has  made  a  calculation  for  a  compact  tension  speci¬ 
men.  His  results  are  essentially  the  same 


Under  a  fatigue  load,  as  the  applied  load  Increases  from  K  fche  stresses 
at  any  point  within  the  cyclic  plastic  zone,  r  ^  increase.  The  plastic  defor¬ 
mation  at  the  point  takes  place  when  the  effective  stress  increases  by  the 
amount  of  2a '  .  ,  instead  of  aY  as  in  the  static  case.  Therefore,  one  can  justi¬ 

fiably  use  Equation  (11)  to  calculate  fatigue  crack  growth  rate  with  the  substi¬ 
tutions  of  da/dN,  AK,  and  2  ct ,  .  for  Aa  ,  K,  and  a  respectively.  After  these 
substitutions,  one  obtains  '  '  uZ 
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0.02  (1-v2)  M 


EO, 


Y(C) 


0.02 


AJ 

°Y(C) 


0.018 


Ea, 


Y(C) 


(12) 


Bates  and  Clark  [15]  related  fatigue  striation  spacing  with  AK  for  a  number 
of  materials.  Taking  striation  spacing  as  da/dN,  we  have 


da 

dN 


6(fV 


(13) 


From  the  available  literature,  Hahn  et  al.  [16]  have  correlated  da/dN  with 
AK  for  a  number  of  steels.  They  found  that 


da 

dN 


(14) 


Pook  and  Frost  [17]  have  found  a  similar  relation 


Barsom  [18]  tested  the  fatigue  crack  growth  of  a  number  of  steels.  He  found 


—  ■  0.66  x  10~®  AK^’^  (inches/cycle) 


(16) 


According  to  Bates  and  Clark,  Hahn  et  al.,  and  Pook  and  Frost,  da/dH  can  be 
written  as 


da 

dN 


C(f, 


(17) 


With  °y(C)^  =  and  E  =  ^07  x  MPa  “  30  x  10^  psi  for  steel,  the  empiri¬ 
cal  equations  and  the  theoretical  equation  are  plotted  in  Figure  15.  The  agree¬ 
ment  is  very  good. 
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Fig.  15.  Comparison  of  the  empirical  fatigue  crack  growth 
rates  in  steels  with  that  predicted  by  the  un¬ 
zipping  model.  Ref.  [ll] 


The  measured  crack  growth  data  of  various  steels  do  not  seem  to  vary  much. 
It  is  known  that  the  cyclic  yield  strengths  of  steels  do  not  vary  much  because 
the  annealed  soft  steels  usually  cyclically  harden,  and  hardened  steels  usually 
cyclically  soften.  The  value  of  the  ratio  <Jy^/E  does  not  vary  much  from  one 


steel  to  another.  If  their  crack  growth  is  controlled  by  the  deformation  pro¬ 
cess,  their  growth  rates  should  not  vary  much  either. 

Kobayashi  et  al.  [19]  measured  the  striation  spacings  of  a  variety  of 
metals  and  alloys.  The  measured  striation  spacings  are  related  to  AJ/a^Qw  as 

shown  in  Figure  16.  For  a  number  of  materials,  ^f^ow  <  ^  0Y(C)  <  ^  ^flow'  Two 


Fig.  16.  Comparison  of  striation  spacing  with  that  prediced 
by  the  unzipping  model.  Ref.  [20] 


lines  in  Figure  16  are  calculated  according  to  Equation  (12)  with  2  "  °f low 

for  one  line  and  °Y(c)  =  °flow  for  the  other*  Aga^n  tfie  agreement  is  very  good. 

More  detailed  discussions  of  the  modelling  and  the  comparison  with  the  experi¬ 
mental  results  are  given  by  Kuo  and  Liu  [9],  Liu,  Yang,  and  Kuo  [3],  and  Liu  and 
Kobayashi  [20]. 

In  the  unzipping  model,  the  crack  growth  rate  is  much  less  than  that  given 
by  the  Dugdale  strip  yielding  model.  As  a  crack  tip  is  deformed  and  blunted 
during  a  rising  load,  the  crack  tip  opens  up  and  the  tip  grows  forward.  However, 
only  part  of  the  blunting  crack  opening  displacement  contributes  to  crack  growth. 

Figure  17a  shows  slip  band  A  at  a  crack  tip.  As  K  increases,  the  shear  de¬ 
cohesion  at  A  causes  crack  tip  blunting  and  crack  growth  as  shown  in  Figure  17b. 
As  the  crack  tip  grows  during  a  rising  K,  the  slip  band  A  is  left  behind.  The 
slip  band  A  will  continue  to  be  active.  Once  behind  the  crack  tip,  the  shear 
decohesion  at  A  will  blunt  the  crack  tip  but  will  not  contribute  toward  crack 
growth  as  shown  in  Figure  17c. 

Let  6KX  -  (1^  -  0);  6K2  -  (Kj  -  Kj);  . . . . ;  6Kr  -  <Kn  -  K^).  Along  the 
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Fig.  17.  Shear  decohesion,  crack  tip  blunting,  and  crack 
growth . 

crack  path,  n  slip  bands  will  participate  in  the  unzipping  process.  After  6Kp, 
the  crack  tip  reaches  the  slip  band  (P  +1).  During  each  of  the  SKjj  increments, 
shear  decohesions  take  place  on  all  of  those  slip  bands  at  or  behind  the  crack 
tip.  All  of  these  shear  displacements  will  blunt  the  crack  tip,  but  only  the 
opening  displacement  on  the  slip  band  at  the  crack  tip  contributes  to  crack 
growth.  Let  (SCTOD^  be  the  crack  tip  opening  displacement,  which  takes  place  on 

slip  band  m  during  the  increment  6K  .  6CT0D  are  tabulated  below. 

P  mp 
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.  .  .  Band  n 

<5^ 

6CT0D1]L 
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6CTOD2n<  .  .  . 
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nn 

The  crack  growth  is  the  sum  of  the  diagonal  terms.  The  crack  blunting  and 
Dugdale  CTOD  include  not  only  the  off-diagonal  terms  but  also  the  shear  deforma¬ 
tion  along  the  crack  flank  which  is  behind  the  original  crack  tip  at  the  be¬ 
ginning  of  a  loading  cycle.  Therefore,  the  crack  growth  rate  calculated  with 
the  Dugdale  model  is  nearly  ten  times  greater  than  the  experimental  values. 

In  the  calculation,  the  slip  step  size,  6Aa  ,  is  assumed  to  be  very  small. 

uz 

The  model  is  valid  even  if  the  slip  step  is  equal  to  one  atomic  spacing. 


A  MODEL  FOR  CRACK  GROWTH  THRESHOLD,  AKth 

Slip  bands  can  easily  be  seen  on  the  surface  of  a  plastically  deformed  crys¬ 
tal.  Intense  plastic  deformation  takes  place  on  certain  parallel  slip  planes, 
while  the  layers  of  the  crystal  between  these  active  slip  planes  remain  practi¬ 
cally  undeformed,  as  illustrated  in  Figure  J2.  In  aluminum,  the  thickness  of 
the  undeformed  crystal  layer,  is  about  200  A  [l].  Let  the  minimum  thickness  of 

the  crystal  layer  be  Z  .  When  the  unzipping  crack  increment  Aa  is  less  than 

m  uz 

Z  ,  the  applied  load  will  reverse  the  slip  direction  before  the  crack  front 
m 

reaches  the  intersecting  conjugate  slip  plane  8  as  shown  by  the  dashed  line, 
position  2  in  Figure  18b.  Upon  unloading,  however,  the  crack  tip 


Fig.  18.  Cyclic  creep  shear  decohesion  when  Aa^z  <  ^  .  Successive  positions 

of  the  upper  crack  surface  during  a  stress  cycle  are  Indicated. 

The  additional  shear  decohesion  between  positions  2  and  4  is  caused 
by  cyclic  creep. 


will  not  return  to  its  original  location.  The  mean  shear  stress  on  the  slip 

plane,  (x  -  X  .  )/2  lies  in  the  direction  to  open  the  crack  tip.  During  the 
max  min 

next  loading  cycle,  the  mean  shear  stress  will  cause  "cyclic  creep"  slip, moving 
the  crack  tip  a  little  closer  to  the  conjugate  slip  plane  8,  position  4  in 
Figure  18b.  Additional  shear  decohesion  between  positions  2  and  4  is  caused  by 
cyclic  creep.  This  cyclic  creep  slip  motion  will  continue  during  the  subsequent 
loading  cycles,  until  the  crack  front  reaches  the  slip  plane  8.  Then,  the  un¬ 
zipping  process  will  switch  to  the  conjugate  plane.  A  large  number  of  cycles  is 
needed  to  propagate  the  crack  by  the  amount  of  Z^.  Therefore,  when  crack  growth 
rate  is  less  than  Z^,  the  growth  rate  decreases  drastically  with  AK.  This  is 
known  commonly.  This  crack  growth  rate  transition  takes  place  at  da/dN  £  Z  and 
at  AK  -  AKt.  “ 

Below  the  transition  point,  the  size  of  the  crack  growth  "step"  remains  con¬ 
stant,  equal  to  Z^.  But  it  takes  a  large  number  of  cycles  to  propagate  a  crack 
through  Iq.  Hertzberg  and  Mills  [2l]  have  observed  constant  crack  growth  step 
in  aluminum,  stainless  steel,  brass,  and  Ni-Cb  eutectic  composite.  Constant 
growth  step  size  were  observed  over  two  decade  range  of  crack  growth  rates. 
Similar  observations  have  been  made  by  Gell  and  Leverant  [22] . 


The  sharp  decrease  in  crack  growth  rate  as  AK  decreases,  often  occur  in  the 
neighborhood  of  200  A/ cycle,  10”6  inches/cycle),  a  value  close  to  the  observed 

Ig  for  aluminum.  The  crack  growth  step  size  observed  by  Hertzberg  and  Mills  is 
likely  to  be  2lg  consisting  of  one  step  on  one  slip  plane  and  another  on  the  con¬ 
jugate  slip  plane. 


The  transition  point,  AKt,  is  given  by 


AK‘ 
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0.018 


EO 


y(c) 


(18) 


When  cyclic  creep  Cakes  place,  the  mean  stress  Is  relaxed.  Whe  the  mean 

stress  is  relaxed,  two  effects  will  increase  shear  decohesion.  First,  as  the 

mean  stress  relaxes  to  zero,  a  oart  of  the  elastic  strain  at  K  is  convert- 

max 

ed  to  plastic  strain  and  shear  decohesion.  Second,  when  the  stress  relaxation 

takes  place,  the  flow  stress  at  K  is  reduced.  This  reduction  in  flow  stress 

max 

will  increase  the  plastic  zone  size  and  increase  plastic  deformation  and  shear 

decohesion.  A  rigorous  evaluation  of  the  stress  relaxation  and  cyclic  creep 

shear  decohesion  is  difficult  especially  if  the  discrete  deformation  mechanism 

is  taken  into  consideration.  However,  the  end  results  of  these  two  effects  can 

be  considered  as  the  reduction  of  c,F/„.  in  the  calculation  of  Aa  . 

Y(C)  uz 

When  stress  relaxation  takes  place,  the  flowstress  at  K  decreases.  Fig- 

max 

ure  19  illustrates  the  stress  relaxation  at  a  point  cycled  within  a  fixed  strain 


*  Fig.  19.  Relaxation  of  mean  stress  and  stablized  hysteresis 

loop.  Successive  stress  reversals  are  numbered. 

range,  Ae.  The  omax  decreases  cycle  by  cycle  until  the  mean  stress  is  relaxed  to 
zero  and  is  nearly  the  half  of  the  value  at  the  beginning  of  the  cycle  load¬ 

ing.  Figure  19  illustrates  the  stress  relaxation  where  plastic  deformation  takes 
place  during  both  loading  and  unloading  half  cycles.  In  this  case,  the  rate  of 
relaxation  is  rather  fast.  However,  the  stress  relaxation  within  the  cyclic 
plastic  zone  is  coupled  to  the  slow  relaxation  within  the  monotonic  plastic  zone, 
where  relaxation  takes  place  only  when  the  stress  deviates  slightly  from  the  elas¬ 
tic  deformation  during  the  re-loading  half  cycle.  The  rate  of  unzipping  cyclic 
creep  shear  decohesion  is  directly  related  to  the  overall  relaxation  rate.  When 
the  mean  stress  is  relaxed  to  zero,  the  flow  stress  is  reduced  approximately  by 
half.  Thus  we  have 


Aauz  •  °'018  |(o. 
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When  the  cyclic  creep  shear  decohesion  is  taken  into  consideration,  if  Aa  is 

uc 


still  less  than  the  crack  tip  will  move  back  and  forth  between  two  neighbor¬ 
ing  conjugate  slip  bands,  and  the  crack  will  remain  dormant.  Therefore  AK  be¬ 
comes  AK^.  Combining  Equations  (18)  and  (19) ,  we  have 


0.7  AK 


AJ  =  0.5  AJ_ 
th  t 


Figure  20  shows  the  data  of  A533B  steel  [23].  The  values  of  AJ  and  AJ 

2  2  ^  th 

are  18  in-lb /m  and  36  in-lb /m  ,  and  the  ratio  AJ  /AJ  ,  is  0.5  as  given  in  Equa¬ 
tion  (20).  ,0., - , - r_ - r-TC-~r-l 


OR  (4K(*/E.  IN-LB* 

Fatigue  crack  growth  rate  of  A533B  steel  versus  AJ. 
Data  from  Ref.  [23] 


The  cyclic  creep  shear  decohesion  model  for  AK^  takes  the  deformation  pro¬ 
cess  into  consideration.  It  is  well  known  that  AK  .  nis  environmentally  sensi¬ 
tive.  The  chemical  effects  must  be  superimposed  onto  the  deformation  mechan- 


The  details  remain  to  be  studied.  However,  the  relation  AK 


0.7  AK 


-  — - - —  — - —  — — ,  - - “  t 

seems  valid  for  a  number  of  materials  tested  in  detrimental  chemical  environ¬ 
ments. 


THE  GROWTH  OF  FATIGUE  MICRO-CRACKS 

When  a  crack  is  very  small,  the  applied  stress  is  often  high.  It  is  often 
close  to  or  even  higher  than  the  yield  stress  of  the  material,  in  order  to  keep 
AK  above  AK^.  In  this  situation,  the  linear  elastic  fracture  mechanics  is  no 

longer  applicable.  Yang  and  Liu  [10]  have  used  the  finite  element  method  to 
model  the  unzipping  crack  growth  process  in  a  single  edge-cracked  plate  with  an 
applied  stress  up  to  the  yield  stress.  The  unzipping  crack  increment  is  given  by 


At  various  K- levels,  the  quantity  dAa  /dK  was  calculated,  and  Aa  was  obtained 

uz  US 

by  numerical  Integration.  The  calculated  crack  growth  rate  is 


(23) 
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where  Ao  is  the  applied  stress  range  and  AK^^  is  the  apparent  AK  value  obtained 

from  the  linear  elastic  solutions.  In  the  case  of  small  scale  yielding,  the 
ratio  Aa/aY^  approaches  zero,  and  Equation  (23)  is  reduced  to  Equation  (12). 

AKg££  is  the  AK-value  corresponding  to  Aa^z  in  the  small  scale  yielding  case. 

Combining  Equations  (12  and  23),  we  have 
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(24) 


In  the  intermediate  crack  growth  rate  region,  Aa^z  corresponds  directly  to 

da/dN.  In  the  lower  AK  region  near  the  threshold,  da/dN  is  much  less  than  Aauz, 

because  of  the  discrete  nature  of  the  crack  growth  process.  However,  in  this  re¬ 
gion,  we  can  use  Aa^z  as  a  physical  parameter  that  characterizes  the  crack  tip 

shear  decohesion  process  as  illustrated  in  Figure  15.  This  is  analogous  to  the 
use  of  elastic  AK  to  characterize  crack  tip  stress  and  strains  even  if  plastic 
deformation  takes  place. 


When  a  nucleated  micro-crack  is  at  the  threshold  of  crack  growth,  AK^^  be¬ 
comes  AK  ,  ,  and  the  applied  stress  amplitude  (Ao/2)  becomes  the  endurance  limit, 
th 

°e*  AKth  iS  related  t0  uz (th)  by  E9uation  (12)>  Aauz(th)  ls  used  here  as  a 
physical  parameter  that  characterizes  the  shear  decohesion  process  at  a  crack  tip 

for  crack  growth  threshold.  With  AK  =1.1  Aa/rra  for  a  single-edge-cracked 

appa 

plate  and  "a"  as  the  crack  nucleus  size,  an,  we  have 
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Equation  (25)  relates  endurance  limit  to  crack  nucleus  size,  a^  and  Aauz^th^ 

The  crack  nucleus  size  is  dictated  by  the  nucleation  process.  For  example,  if 
a  fatigue  crack  is  nucleated  by  an  inclusion  particle,  an  must  be  related  to  the 
size  of  the  particle. 

FATIGUE  CRACK  GROWTH  IN  A  FERRITIC-MARTENSITIC  TWO  PHASE  STEEL 


The  fatigue  crack  growth  behaviours  of  martensitic-ferritic  steel  have  been 
studied  [24,25,26,27],  Fatigue  cracks  in  such  steels  are  often  nucleated  in  the 
soft  ferrite  region.  If  the  applied  stress  is  high  enough,  the  nucleated  cracks 
will  propagate  across  the  martensite  to  the  final  failure. 

Yang  and  Liu  [ll]  analyzed  the  unzipping  crack  growth  process  by  using  the 

finite  element  method.  The  two  phase  structure  is  simulated  by  a  thin  layer 

of  ferrite  over  a  martensite  core  as  shown  in  Figure  21.  The  quantity,  dAa  /dK 

uz 

is  calculated  at  various  K-levels,  and  the  crack  growth  rate  ls  obtained  by  nu¬ 
merical  integration.  The  results  of  the  calculation  are  shown  in  Figure  (22) 
for  three  stress  ranges:  A 0  ■  690,  550,  and  410  MPa  (100,  80,  and  60  ksi).  The 
unzipping  crack  increment,  Aauz  can  be  considered  as  the  crack  growth  rate.  It 

can  also  be  used  as  a  physical  parameter  to  characterize  crack  tip  shear  deco¬ 
hesion  process. 
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Fig.  21.  Finite  element  mesh  represents  a  two-phase 
martensitic-ferritic  steel. 
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Fig.  22.  Unzipping  increment  obtained  by  FQi  simulation,  [ll] 


At  a  given  applied  Ao-level,  before  the  plastic  zone  touches  the  ferrite- 

martensite  interface  boundary,  Aa  increases  with  crack  length  a.  When  the 

uz 

plastic  zone  reaches  the  interface,  the  hard  martensite  constrains  the  plastic 
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deformation  in  the  ferrite.  Therefore,  the  rate  of  increase  of  Aauz  decreases. 
A  maxima  is  reached,  and  Aauz  decreases  as  the  crack  tip  approaches  the  inter¬ 
face.  This  variation  of  the  calculated  Aauz  agrees  with  the  crack  growth  data 
in  Figure  23.  The  growth  rate  increases  initially.  However,  the  growth  rate 


NUMBER  OF  CYCLES. N 

Fig.  23.  Fatigue  crack  growth  in  a  ferritic-martensitic  steel. 

Ref.  [23] 

decreases  as  the  crack  tip  approaches  the  interface.  If  the  applied  stress  is 
low,  and  AK^^  is  below  AK^,  a  crack  may  even  be  stopped  within  the  ferrite  re¬ 
gion.  AK^  is  directly  related  to  Therefore,  the  threshold  condition 

can  also  be  characterized  by  ^auz(t^)-  T^e  qualitative  agreement  between  the 

calculation  and  the  experiment  is  very  good.  An  improved  calculation,  using  the 
actual  material  properties  of  the  ferrite  and  the  martensite,  is  needed  for 
quantitative  comparisons. 

A  QUANTITATIVE  ANALYSIS  OF  MICRO-STRUCTURE  AND  FATIGUE  LIMIT 

In  order  to  have  fatigue  failure  in  a  two  phase  steel,  a  crack  has  to  be 
nut-’eated,  and  the  applied  stress  has  to  exceed  the  threshold  for  the  propaga¬ 
tion  of  the  nucleated  micro-crack.  The  crack  nucleus  in  the  soft  ferrite  has  to 
overcome  the  constraint  of  the  hard  martensite;  the  crack  has  to  penetrate  the 
ferrite-martensite  interface.  Once  in  the  martensite,  the  crack  has  to  overcome 
the  cracking  resistance,  AK^,  of  the  martensite.  The  fatigue  limit  of  the  steel 

can  be  controlled  by  any  one  of  these  processes. 

If  the  crack  nucleus  sizes,  the  crack  nucleation  stresses  the  thresholds 
for  fatigue  crack  growth  in  ferrite  and  martensite  and  the  cyclic  yield  strengths 
of  ferrite  and  martensite,  and  in  addition  the  thickness  of  the  ferrite  layer 
are  known.  The  fatigue  limit  of  a  two  phase  steel  can  be  evaluated  quantita¬ 
tively. 

First,  consider  the  ferrite  alone  without  the  complication  of  the  marten¬ 
site.  Let  the  stress  for  crack  nucleation  be  on.  The  value  of  on  and  will 
give  us  the  value  of  AK^^.  Equation  (24).  If  AK^^  >  (ferrite)  t*'e  crack 

will  continue  to  grow  until  failure.  Therefore,  the  fatigue  limit  oe  is  equal 
to  on,  and  oe  is  controlled  by  the  crack  nucleation  process. 


If  AK  r,  at  a  is  less  than  AlC  ,  ,,  _  s,  the  nucleated  crack  will  stay 

eff  n  th(ferrite)  3 

dormant,  until  the  applied  stress  is  increased.  Therefore,  o  >  o  and  the  fati¬ 
gue  limit  is  controlled  by  the  crack  propagation  process. 


The  relation  between  crack  nucleus  size  and  fatigue  limit  is  given  by  Equa¬ 
tion  (25).  For  the  ferrite,  we  assume  =  380  MPa,  =  3.8  x  10"  , 


and  Aa  .  .  .=  28  x  10-6  mm,  which  corresponds  to  AK  ,  =  11  MPav^n.  The  calcu- 
uz(th)  th 

lated  relation  between  an  and  ae  is  plotted  in  Figure  24.  ae  changes  from  200 

to  300  MPa  as  a  varies  from  0.085  to  0.03  mm. 


n 


The  curve  for  the  martensite  is  also  plotted  in  the  figure.  For  the  mar¬ 
tensite,  we  assume  =  23  x  10~6  mm,  which  corresponds  to  AK^  = 

16.5  MPa»4^  .  ay(c)  =  *®35  MPa>  anc3  E  =  207,000  MPa.  The  fatigue  limit  of  the 

martensite  changes  from  600  to  800  MPa  as  the  crack  nucleus  size  varies  from 
0.02  to  0.01  mm. 


AK 


Once  a  crack  grows  into  the  martensite,  if  the  applied  AK  is  less  than 

the  crack  will  stop  propagating.  The  fatigue  limit  is  thus 


th(martensite) ’ 

controlled  by  the  cracking  resistance  of  the  martensite.  With  the  crack  length 
equal  to  the  ferrite  layer  thickness,  the  fatigue  limit  can  be  calculated  with 
Equation  (12  or  25).  For  a  single  edge-cracked  plate,  and  for  ^^(nartensite)’ 

16.5  MPa  the  fatigue  limit  of  the  two  phase  steel  is  plotted  in  Figure  25. 


The  fatigue  limits  controlled  by  crack  non-propagation  in  both  martensite 
and  ferrite  are  also  plotted  in  the  figure.  For  the  ferrite,  0*  “  300  MPa  and 
an  “  0*03  mm;  and  for  the  martensite,  oe  ■  800  MPa  and  Sq  ■  0.01  mm. 

When  a  crack  tip  grows  close  to  the  ferrite-martensite  interface,  the  un¬ 
zipping  crack  increment,  Aa^  decreases  as  shown  in  Figure  22.  When  the  shear 
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Fig.  25.  Fatigue  limit  and  ferrite  domain  size. 

decohesion  increment  decreases  below  ^auz(th)  £°r  crack,  growth  threshold,  the 

crack  stops  in  the  ferrite  region  as  shown  by  the  three  specimens  at  the  lower 
stress  levels  in  Figure  23. 

The  unzipping  crack  increment  calculation  as  shown  in  Figure  22  could  be 
used  to  construct  the  relation  between  the  fatigue  limit  of  the  steel  and  the 
ferrite  domain  size,  when  martensite  constraining  is  the  controlling  process. 
However,  a  much  finer  finite  element  mesh  is  needed  to  establish  the  quantita¬ 
tive  relation.  The  line  in  Figure  25  is  a  schematic  line. 

Figure  25  gives  an  overall  view  of  the  fatigue  limit  of  the  steel.  The 
effects  of  crack  nucleus  size  in  the  ferrite  and  the  martensite,  the  crack 
growth  thresholds  of  the  ferrite  and  the  martensite,  the  stress  necessary  for 
crack  nucleation,  and  the  ferrite  domain  size  are  all  intimately  related  to  the 
fatigue  limit.  Such  an  analysis  will  help  to  optimize  the  micro-structures  in 
order  to  achieve  the  best  fatigue  strength. 

This  analysis  illustrates  the  need  to  study  the  effects  of  large  scale 
yielding,  the  discrete  nature  of  plastic  deformation,  and  the  material  inhomo¬ 
geneity  in  order  to  establish  the  quantitative  relation  between  fatigue 
strength  and  micro-structure. 

DISCUSSIONS 

Fatigue  crack  growth  rate  are  often  correlated  with  AK.  If  the  specimen  is 
large  enough  and  the  applied  load  is  low  enough  so  that  the  condition  of  small 
scale  yielding  prevails,  AK  is  capable  of  characterizing  crack  tip  stresses  and 
cyclic  plastic  strains.  In  other  words,  a  given  value  of  AK  corresponds  unique¬ 
ly  to  a  state  of  crack  tip  stresses  and  strains.  If  fatigue  crack  growth  is  the 
result  of  cyclic  stresses  and  cyclic  strains  at  the  tip,  da/dN  must  be  related 
uniquely  to  AK.  When  the  condition  of  small  scale  yielding  is  not  met,  a  given 
value  of  AK  may  result  in  a  number  of  different  states  of  crack  tip  stresses  and 
strains.  Therefore,  AK  is  no  longer  able  to  correlate  with  da/dN.  In  the  case 


of  large  scale  yielding,  both  Aauz  and  AJ  are  applicable. 

AK  and  AJ  are  applicable  as  long  as  they  can  characterize  crack  tip  stresses 
and  strains.  Both  AK  and  AJ  are  used  as  indirect  parameters,  which  do  not  infer 
any  crack  growth  mechanisms.  Nor  the  growth  mechanism  need  to  be  known  in  order 
to  make  such  correlations. 

In  the  case  of  the  ferritic-martensitic  two  phase  steel,  neither  AK  nor  AJ 

is  capable  to  characterize  the  crack  tip  stresses  or  strains.  Therefore,  neither 

AK  nor  AJ  can  be  used  to  correlate  crack  growth  rate.  Aa  is  equal  to  da/dN; 

uz 

therefore  it  is  a  direct  parameter  for  crack  growth  correlation.  A  direct  para¬ 
meter,  if  known,  is  always  preferred  because  it  has  a  much  wider  range  of  appli¬ 
cability. 

Mechanical  metallurgy  studies  the  mechanical  properties  as  affected  by 
microstructures.  In  the  stress  analysis,  neither  discrete  plastic  deformation 
nor  material  inhomogeneity  are  taken  into  consideration.  Fatigue  crack  growth 
is  the  direct  result  of  the  discrete  plastic  deformation  process.  A  quantita¬ 
tive  relation  can  be  derived  directly  only  if  the  discrete  deformation  process 
is  taken  into  consideration. 

It  is  suggested  to  develop  a  new  field  of  study  of  metallurgical  mechanics, 
which  takes  inhomogeneity,  anisotropy,  and  the  discrete  deformation  process  in¬ 
to  account  in  the  stress  analysis.  Then,  use  the  calculated  local  stresses  and 
strains  to  correlate  with  microstructures.  This  paper  illustrates  such  studies 
on  the  effects  of  discrete  slip  process  and  the  effects  of  material  inhomogeniety 
on  crack  growth.  The  recent  developments  in  mechanics  and  metallurgical  studies 
set  the  stage  for  fruitful  efforts  and  rapid  progress  in  this  area. 

CONCLUSIONS 

1.  Cyclic  plastic  deformation  at  a  crack  tip  may  cause  a  fatigue  crack  to  grow. 
In  the  case  of  small  scale  yielding,  AK  is  able  to  characterize  the  cyclic  plas¬ 
tic  deformation  at  a  crack  tip,  and  AK  can  be  used  to  correlate  with  da/dN  suc¬ 
cessfully. 

2.  The  unzipping  model  of  fatigue  crack  growth  agrees  with  the  measured  crack 
growth  rate  and  striation  spacing  measurements. 

3.  The  unzipping  model  is  extended  to  micro-cracks  and  the  ferritic-marten¬ 
sitic  two-phase  steel. 

4.  AK  and  AJ  are  indirect  parameters  to  correlate  da/dN. 

5.  Neither  AK  nor  AJ  can  be  used  to  study  the  micro-structural  effects  on 
da/dN.  In  this  case,  a  direct  parameter,  such  as  Aa^  is  needed. 

6.  A  model  of  crack  growth  threshold  is  proposed.  The  model  predicts  cor¬ 
rectly  the  crack  growth  behavior  in  the  near  threshold  region. 

7.  A  quantitative  analysis  on  raicrostructure  and  fatigue  limit  is  illustrated. 
Such  quantitative  analysis  will  help  to  optimize  the  micro-structures  to  achieve 
high  fatigue  strength. 
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